


STRUCTURAL DYNAMIC ANALYSIS 01 
ROTOR BLADES WITH SWEPT TIPS 


M. VENU GOPAL 




INEERIIMOi 





STRUCTURAL DYNAMIC ANALYSIS OF 
ROTOR BLADES WITH SWEPT TIPS 


A Thesis Submitted 
in Partial Fulfilment of the Requirements 
for the Degree of 
MASTER 0E TECHNOLOGY 


by 

M.VENU GOPAL 



DEPARTMENT OF AEROSPACE ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 


MAY, 1999 



CERTIFICATE 



It is certified that the work contained in the thesis entitled ’’STRUCTURAL 
DYNAMIC ANALYSIS OF ROTOR BLADE WITH SWEPT TIPS ”, by 
M.VENU GOPAL has been carried out under my supervision and that this work 
has not been submitted elsewhere for a degree. 



Dr. C.VENKATESAN 
Associate Professor 
Dept, of Aerospace Engineering 
I.I.T. Kanpur 


May, 1999 



o 1 JUH 1999 , 

CENTRAL LifcSKAfO 

$ J. T., KANPUR 


<•«** 12MU 


ACKN OWLEDGEMENTS 


I take this opportunity to thank my thesis supervisor Dr. C.Veukatesan . for 
h,s inspiration and valuable suggestions in carrying out this work. I consider myself 

fortunate to have worked under him for the frppdnm a a 

eedom and constant encouragement he 

gave throughout this work. 

I express my special thanks to my friends Voi;o u ^ ,• 

* as VallBa bu, Kahprasadanna, Ravianna, 

Ramakoti, Vikky, Vinay, Narsi, Lala, Santhosh Ru;„ r>^ „ . 

’ osl1 ’ Wb Bhargav, Gopal, Janardhan, 

Sharath, Giri, Chandv 2 , Chadra sekhar, Kkis q r r> • 

kis, b.B.Raoanna, Rajananna,Prasanth 

for their tremendous help in bringing this work to the present form. 

I would like to thank Mr. Rameshwar, Siva Anil OcJ a 

j Anil, Sada and many other friends 

for their suggestions during my course of study at IIT K 

I have no words to express my thanks to mv narpty+c „ t -i 

y parents and my family members who 

have been constant source of inspiration to me T w ioii +u 1 n c ■ , 

me - 1 wish to thank all my friends and 

well washers who made mv stay at IIT KANPTIP ,, 

J !Nru ^i memorable and pleasant. 


M.VENU GOPAL 



Contents 


Abstract iv 

List of Figures v 

List of Tables vii 

List of Symbols viii 

1 INTRODUCTION 1 

1.1 STRUCTURAL MODELLING OF ROTOR 

BLADE 2 

1.2 OBJECTIVES 4 

2 ROTOR BLADE MODEL AND ASSUMPTIONS 5 

2.1 BASIC ASSUMPTIONS 5 

2.2 ORDERING SCHEME 6 

3 COORDINATE SYSTEMS 8 

3.1 HUB FIXED INERTIAL SYSTEM-R 8 

3.2 HUB FIXED MOVING SYSTEM-H 9 

3.3 HUB FIXED ROTATING SYSTEM-1 9 

3.4 ROTATING SYSTEM-2K 10 

3.5 PRECONED, ROTATING SYSTEM-3K 10 


i 



11 


3.6 

PREDROOPED, PRESWEPT, PITCHED, 



BLADE-FIXED ROTATING SYSTEM 4K 

. ... 10 

3.7 

UNDEFORMED ELEMENT COORDINATE SYSTEM-e . . . 

. . . . 11 

3.8 

ROTATING, BLADE-FIXED SYTEM .... 

. . . . 12 

3.9 

COORDINATE SYSTEM-6K - - • 

. . . . 13 


4 KINEMATICS 

4.1 POSITION VECTOR OF A POINT 

4.2 ANGULAR VELOCITY VECTOR 

4.3 VELOCITY AT POINT ’P’ . • ■ ‘ 


5 EQUATIONS OF MOTION FOR ROTOR BLADE 

5.1 KINETIC ENERGY OF BLADE • 

5.2 STRAIN ENERGY OF BLADE • • 

5.2.1 Strain Energy • • ' 

5.2.2 Explicit Strain-Displacement Relations 

5.2.3 Stress-Strain Relations • • 

5.2.4 Strain Energy Variation • 


• FORMULATION OF ELEMENT MATRICES ASSOCIATED WITH 
KINETIC AND STRAIN ENERGY VARIATION 34 

6 1 FINITE ELEMENT DISCRETIZATION 34 

w(rpmnrc ASSOCIATED WITH KINETIC ENERGY 

6.2 ELEMENT MATRICES ASSUW* 

37 

VARIATION • 

gg 

6.2.1 Mass Matrix [M) 14xl4 • • ' ' 

41 

6.2.2 Matrix [M c ]i4xi4 • • • * ' ‘ ’ 

6.2.3 Centrifugal Stiffness 44 

47 

6.2.4 Vector [V L ]i4xi 

. 48 

6.2.5 Matrix [M Ji4x3 • • * ' ' 



Ill 


6.2.6 Matrix [A/ 2 ]i 4x3 50 

6.2.7 Matrix [M 3 ] 14x3 52 

6.2.8 Matrix [M 4 ]i 4x3 56 

6.2.9 Vector [V 14xl ] 58 

6.2.10 Nonlinear Vector {V r,fVi ]i 4xl 59 

6.3 ELEMENT MATRICES ASSOCIATED WITH STRAIN ENERGY 

VARIATION 60 

6.3.1 Linear Stiffness Matrix [ K E ] 60 

6.3.2 Simplification of Linear Stiffness Matrix [ K E ] 62 

6.3.3 Nonlinear Stiffness Vector [ F E ] 67 

6.3.4 Treatment of Nonlinear Terms Associated with Axial Strain at 

Elastic Axis 68 

7 RESULTS AND DISCUSSION 72 

7.1 VALIDATION 73 

7.2 EFFECTS OF TIP SWEEP AND ANHEDRAL 74 

8 CONCLUDING REMARKS 76 

References 107 

A Local to global coordinate transformation for tip element 111 

B Natural frequencies and coupling measures of rotor blade with swept 
tips 


113 



IV 


ABSTRACT 

This thesis presents the formulation of the structural dynamic equations of motion of 
a general rotor blade configuration with swept tip. The equations of motion have been 
derived using Hamilton’s principle and Rayleigh Ritz finite element discretisation of 
the beam. First the formulation is validated by comparing the results of the present 
analysis for a straight blade with those available in the literature. Then a systematic 
structural dynamic analysis has been performed to identify the effect of tip sweep 
angles on the natural frequencies and mode shapes of the rotating blade. 

Quantitative measures of lag-torsion and flap-torsion couplings have been defined 
based on the contribution of torsional deformation in the in-plane bending and out-of- 
plane bending modes. The effect of tip sweep on these couplings have been identified. 
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Chapter 1 

INTRODUCTION 


It is well known that the time varying loads on the main rotor system contributes 
significantly to the vibration in helicopter. Therefore, the structural dynamic charac- 
teristics of the rotor blade and also the dynamic characteristics of the fuselage have 
a very strong influence on the vibratory levels in helicopter. Any analytical study 
of vibration in the helicopter requires the development of a dynamic model of the 
helicopter system. The major components of this model are: 


• Rotor blade model 

• Fuselage model 

• Rotor-Fuselage interface model 

The formulation of the rotor blade model requires the development of structural, 
inertial and aerodynamic operators associated with the rotor blade motion. The he- 
licopter fuselage model is represented by an idealized structural dynamic model of a 
three-dimensional structure. The rotor-fuselage interface model must represent both 
the geometry of the interface as well as the aerodynamic interaction in an appropriate 
manner. 

This study deals with the formulation and solution of a structural dynamic (structural 
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and inertia operators) model of a rotor blade with swept tip incorporating all complex 
geometric parameters of a realistic rotor system. 


1.1 STRUCTURAL MODELLING OF ROTOR 
BLADE 

Since helicopter rotor blades are long, slender beams undergoing moderate deforma- 
tions, a nonlinear strain-displacement model is used to describe the coupling effects 
between axial, bending and torsion deformations. Generally, the strains are assumed 
to be small in comparison to unity. Such an assumption is consistent with the de- 
sign requirement based on fatigue life consideration which states that the rotor blades 
must be designed to have an operating strain level well below the elastic limit of the 
blade material. 

The first structural dynamic model for a rotor blade undergoing flap-lag-torsional 
deformations was developed by Houbolt and Brooks[l]. This model neglects the non- 
linear coupling effects between bending and torsion, which were shown by later re- 
searchers, to be important for the rotor blade analysis. 

Several moderate deflection beam theories have been developed [Refs. 2-7] for he- 
licopter rotor blade analysis. These beam models were later used to formulate the 
structural, inertial and aerodynamic operators for the aeroelastic stability analysis 
of rotor blades [8-11]. Structural modelling for composite rotor blade has also been 
developed recently by various investigators [11-20]. The earlier models presented in 
Refs. 2-7, were restricted to treatment of isotropic blades. In general, these models 
did not include the effects of cross-sectional transverse shear and warping. On the 
other hand, the composite blade models include these effects. 

In the eighties, efforts have been made to design rotor blades with modified tips to 
improve the rotor performance and loads. Subsequently several studies (Refs. 21-25) 
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were undertaken to address the influence of tip sweep on the aeroelastic stability and 
performance of rotor blades. A schematic of a rotor blade with tip sweep is shown in 
Fig. 1. Due to sweep (angle A s defined in the plane of rotation) and anhedral (angle 
A a defined in out-of-plane of rotation), rotor blades experience structural bending- 
torsion coupling effects. In addition due to high dynamic pressure, swept tips greatly 
influence the aeroelastic stability and response characteristics of the blade. Thus, tip 
sweep and anhedral angles can be effectively used to tailor the aeroelastic character- 
istics of the blades, there by imporve the noise and vibratory loads of the rotor. 
During the derivation of the equation of motion of the rotor blade with swept tips, the 
formulation of inertia and aerodynamic operators generate a large number of higher 
order nonlinear terms. A consistent ordering scheme is used by assigning orders of 
magnitude to all the nondimensional parameters of the problem, in terms of blade 
bending slope (which is assumed to be of the order of e ). A second order approxima- 
tion assumes that the terms of order e 2 are neglected in comparison to terms of order 
1. In this study a third order approximation has been used, i.e., terms of order e 3 are 
neglected compared to terms of order 1. i.e., 

1 + 0(e 3 ) « 1 

In all the studies reported in Refs. 2-25, only certain geometric parameters such as 
precone, hinge offset and swept tips have been included in modelling the rotor blade 
dynamics. In the present formulation, an attempt has been made to include all the 
complex geometic parameters such as torque offset, blade root offset, precone angle, 
predroop angle, presweep angle, pretwist, tip sweep angle and tip anhedral angle, as 
well as hub motion. 

Th is study is an extension of Ref. 26, which considered all geometric complexities of 
the rotor blade and hub motion, but did not consider tip sweep. 
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1.2 OBJECTIVES 

The objectives of the present study are: 


• To develop a most general beam type finite element model for the rotor blade 
with, swept tips including various geometric parameters like root offset, torque 
offset, precone, predroop, presweep, pretwist and hub motion. 

• To validate the model by comparing the results of this study with those available 
in literature. 




Conduct detailed studies on the dynamic behavior of rotor blade with straight 

of 

and swept tips to determine the effects. bending-torsion coupling due to sweep 

A 


and anhedral. 



Chapter 2 


ROTOR BLADE MODEL AND 
ASSUMPTIONS 


Helicopter rotor blades are long slender beams attached to the hub through a complex 
geometrical and mechanical arrangement. The geometrical parameters describing the 
configuration of the rotor blade-hub system is shown in Fig. 1. The parameter ’a ! 
represents the torque offset which is the distance from the centre of rotation (hub 
center) to the reference axis of the blade, ’ex’ and 5 e 2 ’ refer to the blade root offset 
distance from the hub. P p represents the precone angle defining the orientation of 
the blade pitch control axis with respect to the hub plane. f3 s and Pd represent the 
presweep and predroop angles respectively, representing the inclination of the blade 
reference axis with respect to the pitch control axis. The blade consists of a straight 
portion and a swept tip whose orientation relative to the straight portion is described 
by a sweep angle A s and an anhedral angle A a . 

2.1 BASIC ASSUMPTIONS 

In the forumulation of the dynamic model of the rotor blade with swept tip, several 
assumptions have been made which are given below. 

1. The blade is treated as an elastic beam. 
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2. The blade is modeled by beam type finite elements along the elastic axis of the 
blade. 

3. The speed of rotation (f2) of rotor is constant. 

4. The rotor shaft is rigid. 

5. The blade undergoes moderate deformation in flap, lag, torsion and axial modes. 

6. The blade cross section can have a general shape with distinct shear centre, 
aerodynamic centre and centre of mass. 

7. The blade is assumed to be made of isotropic material but can have nonuniform 
properties along the span. 

2.2 ORDERING SCHEME 

In the formulation of the equations of motion of a rotor blade with swept tip under- 
going moderate deformations, a large number of higher order terms are generated. In 
order to identify and eliminate higher order terms in a consistent manner, an ordering 
scheme is employed. This ordering scheme is based on the assumption that the slopes 
of the deformed elastic blade are moderate, and of order e (0.10< e < 0.20). Orders of 
magnitude are then assigned to various non-dimensional physical parameters govern- 
ing the rotor blade dynamic problem, in terms of (. In the derivation of the governing 
equations, it is assumed that terms of order e are neglected with respect to terms of 
order 1, i.e., 

0(l) + 0(e 3 )«0(l) 

The order of magnitude of various non-dimensional parameters governing this prob- 
lem are given below: 

ORDER 1 


cos<j> k , sin<j>k, A a , A, = 0(1) 
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f = 0( 1 ) 

*£() = £0 = 0(0 
<&() = *0 = 0 ( 1 ) 

ORDER es 
6 G k = O(e^) 


ORDER e 

a £l 22 2 Jit Hi — Off') 

/> f ’ z > c > / > z — 
v' k ,w' k ,(f>,pp,p d ,p s = 0(e) 


ORDER ef 

3 

/m w ,/m cc = 0(es) 

^,^A, 0 y , 0 *=O(e f ) 

ORDER e 2 

Y,u k ,mrjm,m^ m = 0(e 2 ) 


It is important to note that ordering schemes are based on physical understanding 
of the behaviour of actual blade configurations. Hence care must be excercised in 
deleting higher order terms, based on this ordering scheme. 



Chapter 3 

COORDINATE SYSTEMS 


The description of the complex deformation of a rotor blade requires several coordinate 
systems. The transformation relation between quantities referred in various inertial, 
non-inertial coordinate systems is to be established before deriving the equations of 
motion of the rotor blade. The relation between two orthogonal systems Yb Z, and 
Xj, Yj, Zj with e xi , e yi , e zi and e X j, e' y j, e z j as unit vectors along the respective axes can 
be written as 


' e xi ] 


f e x j 

* ^yi 

> = [T u ] - 

| Cyj 

. £zi j 


[ J 


(3.1) 


where the transformation matrix [Ty] can be obtained using the Euler angles 


required to rotate the j-system so as to make it parallel to i-system. The coordinate 
systems used in deriving the equation of motion for the rotor model are described 


below: 


1 HUB FIXED INERTIAL SYSTEM-R 

■ pe coordinate system ’R’, shown in Fig. 2, has its origin at the centre Or of the 
□■perturbed hub. The xr axis is pointing towards the helicopter tail and Zr is 
jointing upwards. The unit vectors along the three axes are e X R,e y R,e Z R. 
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3.2 HUB FIXED MOVING SYSTEM-H 

The coordinate system H shown in Fig. 3, is a body fixed system with its origin 
fixed at the centre of rotor hub 0 H . Prior to perturbational motion of the hub, H- 
system coincides with R-system. If 8 X , 9 y and 9 Z represent the sequential yaw-pitch- roll 
rotations, then the transformation matrix [Thr] can be written as- 

10 0 cos 6 y 0 — sinflj, *1 f cos 9 Z sin 9 Z 0 ' 

[Thr] = 0 cos 9 X sin9 x 0 1 0 -sin 9 Z cos 9 Z 0 (3.2) 

0 — sin0 z cos 9 X _ _ sin 9 y 0 cos 9 y 0 0 1 

Since 9 X , 9 y and 6 Z are assumed to be of order e = , sine and cosine functions can be 
approximated as 

V 

sin 9 as 9, cos 9 « 1 

Substituting the approximation, [Thr] can be simplified as 

1 9 Z -9 y 1 

[Thr] = 9 y 9 z — 9 Z 1 9 X (3.3) 

J x 9 z + e y e y e z -e x 1 

3.3 HUB FIXED ROTATING SYSTEM-1 

The coordinate system 1 shown in Fig. 4, rotates about z\h axis with constant angular 
speed fl of the rotor. Its origin is fixed at hub centre Oh- This system can be obtained 
by rotating H system by an azimuthal angle 'ipk of the k th blade about zm axis. The 
transformation matrix is given as - 


[Tjh] = 


cos tpk sin ipk 0 
— sin ipk cos ‘ipk 0 

0 0 1 


(3.4) 
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where, xpk is azimuth angle of k th blade 


Vk = + (k - 1) — 


and xp = Qt 


3.4 ROTATING SYSTEM-2K 


The coordinate system 2K shown in Fig. 5 is a blade fixed coordinate system which 
rotates with k th blade. The origin of the 2K system is at the location of the blade 
root A (Fig. 1), which is at a distance ae y i + e\e x \ from the hub centre. 1-system and 
2K-system are parallel. • 


3.5 PRECONED, ROTATING SYSTEM-3K 

The system 3K shown in Fig. 6a and Fig. 6b also rotates with blade. This system 
is obtained by rotating 2K-sytem by an angle- /3 P (precone angle) about y 2 k axis. The 
transformation matrix between 2K and 3K system is given as- 


3.6 


pr«] = 


i o p p - 

o 1 o 


(3.5) 


i-Pp 0 1 J 

PREDROOPED, PRESWEPT, PITCHED, 
BLADE-FIXED ROTATING SYSTEM-4K 


The 4K system shown in Fig. 6a and Fig. 6b, is blade fixed system with its origin at 
pitch bearing B of the blade. 

It may be noted that the pitch axis of the blade is along direction and the blade 
reference elastic axis is along the e l4fc direction. While changing the control pitch 
input of the blade, the elastic axis will move on the surface of a cone whose vertex 
angle is described by the angles and as shown in Fig. 6c. 



11 


The 4k system is obtained by the following steps- 

Step-1 Translating the origin of 3K system by a distance ’e 2 ’ along e l3 * 

Step-2 Then rotating the system by an angle -/3 S (preweep angle) about z 3fc axis. 
Step-3 Then rotating the system by an angle -/?d(predroop angle) about y 3k . 

Step-4 Then rotating the system by an angle Oj (pitch input) about xzk- 
The transformation matrix is given as- 

1 — (P s cos 6 1 -I- Pd sin 9 1 ) (Pd cos Oj — P s sin Qi ) 

(P s cos 6j + PdSindi) cost?/ sin 6j 

(— Pd cos Oi + P s sin 6 1 ) —sin 0/ cos Oj 

(3.6) 

3.7 UNDEFORMED ELEMENT COORDINATE 
SYSTEM-e 

The e system, shown in Fig. 7, has its origin at the inboard node of the finite element. 
The vector e xe , is aligned with the beam element elastic axis; while the vectors e ye 
and e ze are cross sectional coordinate axes. For the straight portion of the blade , 
the ( e xe ,e ye ,e ze ) system has the same orientation as (e X 4 k,e yi k,e Z 4 k ) system. For the 
swept-tip element, the e system is oriented by rotating the 4k system about e z4k by 
sweep angle A s , and then about e y4k by the anhedral angle A a . 

The transformation matrix between 4k and e system is given as- 



[ e xe ' 


&x4k 

1 2 

VC 1 

► = [Te4] < 

| &y4k > 

l e 2e J 


[ &z4k J 


For the element in the straight portion of the blade 


f 1 0 0 


[T«4] = 


0 1 0 
0 0 1 


(3.7) 


(3.8) 
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For the swept-tip element 


[t.4] = 


cos A 3 cos A a — sin A s cos A s sin A a 
sin A s cos A a cos A s sin A 4 sin A a 

— sin A a 0 cos A a 


(3.9) 


where, A s is the blade tip sweep angle, positive for backward sweep, and A a is the 
blade tip anhedral angle, positive upward. 


3.8 ROTATING, BLADE-FIXED SYTEM-5K 

The 5k system, shown in Fig. 8 is a cross-sectional coordinate system of the k th blade. 
In the undeformed state of the blade, both e and 5k systems are parallel. But the 
origin of the 5k system is at a distance x k from the origin of the e system. During 
elastic deformation of the blade, the 5k system translates and rotates with the cross- 
section. After deformation, the origin of the 5k system, from the origin of 4K system, 
is at the location given by 

n— 1 

( ^ y {%k *b 'Uk)&xe ~b Vk&ye "b ‘^k^ze (3.10) 

t=l 

The transformation matrix between e and 5k sytem is obtained following a flap -lag 
sequence of rotation. The Euler angles are respectively -P k and corresponding to 
the local slope of the deformed blade in flap and lag directions. The transformation 
matrix is given by: 



COSOc 

sin^ 

0 ' 


cosfi k 

0 

22. 

i 

[T5.) = 

— sin C* 

COS^fc 

0 


0 

1 

0 


0 

0 

1 


. — sin p k 

0 

cos p k _ 


Since the angle £* and P k are of the order 0(e), the transformation matrix can be 
simplified by noting that- 


sin 6 « 6 cosO « 1 
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The Euler angles can be expressed in the terms of the local slope of elastic deforma- 
tion of the blade as 


0k = w'k 
tk = V k 


dw k 

dx 

dvk 

dx 


n— 1 

Note : x = ]T(l e )i + x k 


t=i 


Substituting the above relations in the transformation matrix [T 5e ] yields 


[T 5e ] - 


1 t* w’ k 
-v' k 1 -v' k w' k 
-w' k 0 1 


3.9 COORDINATE SYSTEM- 6 K 


(3.12) 


The 6K system shown in Fig. 9 represents the cross-sectional coordinate system in 
the deformed configuration of the blade. e v — represent the directions of the cross- 
sectional principal axes. 6K-system is obtained by rotating 5K system about e x5k 
through the angle (<j) k + 6c), where do represents the geometric twist angle of the 
cross-section and 4> k represents the elastic twist. 

The transformation relation is given as- 


' e x ' 


'10 o' 


t, 1 

e y 

► = 

0 cos [6 C + 4>k) - sin(0 G + 4> k ) 

< 


UJ 

bk 

_ 0 sin(0 G -1- 4> k ) cos {6 g + <j> k ) 


J 


(3.13) 


Chapter 4 
KINEMATICS 


During operation, the rotor blade undergoes deformation in in-plane bending (lag), 
out-of-plane bending (flap), torsion and axial modes. In addition, the hub centre has 
both translational (R x ,Ry,R z ) and rotational ( 0 x ,6 y ,6 z ) perturbational motion. The 
formulation of inertia operator and aerodynamic operator requires a proper discre- 
tion of kinematics of the blade motion. In this section, an expression for the absolute 
velocity vector at any arbitrary point ’p’ on the blade is derived. 

4.1 POSITION VECTOR OF A POINT 

The position vector of any arbitrary point ’p’ in the n th finite element of the deformed 
blade with respect to the hub centre Oh , is given by 

n-l 

Tp = 0,e y i + €ii x 2k + £2^x3 It + 52(^)t'dx4it + ( x k + Uk)e xe 

i= 1 

+Vke ye + Wkize + t) + C^C (^-1) 


Tr ans for mi ng all the unit vectors to the 4K-system and neglecting the higher order 
terms, the position vector can be written as: 
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n-1 

r p = e xik l[~a{P 3 cos 0/ 4 0 d sin 0,) 4 ei 4 e 2 4 (*«)« 

1=1 

4 cos A s cos A a (:r* 4 u*) 4 t;* sin A s cos A 0 - io* sin A a 
+ 7 ?{~X cos A 4 cos A a cos(0 G 4 0*) - ti4 cos A s cos A a sin(0 G 4 <j> k ) 

- sin A s cos(0 G 4 4>k) 4 cos A s sin A a sin(0 G 4 <j> k )} 

+C{^ cos A s cos A a sin(0 G 4 o k ) - w' k cos A s cos A a cos(0 G 4 <j> k ) 

4 sin A s sin(0 G 4 <j> k ) 4 cos A s sin A a cos(0 G + <j > k )}] 

+e y 4 k l[acos9j + ei(/3 4 cosA; 4 Afsin0 7 - /3 p sin9i) 4e 2 (/? s cos 0/ 

4 A* sin 9j) sin A s (x k 4 u k ) + v k cos A s + r){~v' k sin A s cos A a cos (0 G 4 <p k ) 

—w' k sin A s cos A a sin(0 G 4 <j> k ) 4 cos A s cos(0 G + <j> k ) 

4 sin A s sin A a sin (9 G 4 0*)} + CM sin A s cos A a sin(0 G + <p k ) 

-w' k sin A s cos A a sin (0 G 4 <f> k ) — cos A s si n(9 G 4 <j> k ) 4 sin A s sin A a cos(0 G + 0/t)}] 
+e z4 fc/[— a sin 0/ 4 e\ (-fid cos 0/ + p s sin 0/ - (3 P cos 0/) 4 e 2 (~P d cos 0/ 

+p s sin 0/) 4 cos A s sin A a (x k 4- u k ) 4 v k sin A s sin A a + w k cos A a 
+rj{v' k sin A a cos (0 G + <f> k ) + w' k sin A a sin(0 G 4 <j> k ) 

+ cos A a sin (0 G + (t> k )} + CK sin A 0 cos(0 G 4 (j> k ) - v' k sin A a sin(0 G 4 <p k ) 

4 cos A a cos (0 G 4 <£*)}] 

(4.2) 

In Eq. 4.2 all the length quantities are nondimensionalized with respect to the length 
of the blade 'l'. For the sake of convenience, the nondimensional length quantities are 
referred without an overbar(). The expression given in Eq. 4.2 is general, in the sense 
that for position vector of any point in the straight portion of blade A c and A s are 
set equal to zero. Eq. 4.2 can also be written in symbolic form as ' 
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Tp 1 [t'x&x 4 k ^y^yAk ^z^zAk] 

4.2 ANGULAR VELOCITY VECTOR 

The angular velocity vector u of k th blade consists of three components. They are 


• The constant rotational speed (Q) of the rotor; 

• the rigid body angular velocity (3 r igid) of the hub due to perturbationai rotation 
in roll-pitch-yaw ( 0 x ,6 y and 9 Z ) 

• the angular velocity contribution due to the rate of change of control pitch input 
Q,6j to the blade 

u = £le z u + w T igid + QQi&xZk (4-3) 

where 

frigid = Q@x&xl k "b Q9y&ylk ”b Qflz&zlk (4.4) 

where ( ) indicates differentiation with respect to the nondimensional time ip, 

ip = fit 

Note that 

io = «3ffiO -«£<)- no 

Transforming all the unit vectors of Eq. 4.3 to the 4K-system and neglecting the 
higher order terms, the angular velocity of k th blade can be written in symbolic form 


as: 


LO — \^x^x4k "b LOy^y4k ”b ^z^z4k\ 


(4.5) 
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where 

u) x = [9 X cos ipk 4 6 y sin ip k 4 9j + (P p 4 /3 d cos 0/ 

-&sin0/)] (4.6) 

u> y = [~^* cos Qj sin i>k 4 cos 0/ cos V’jfc + 9 Z sin 9j 
+9 1 (P s cos 9j + Pi sin 9j )] 

<-o z = [9 X sin 9 1 sin ip k - 9 y sin 9j cos ip k 4 9 Z cos 0/ 

4cos0/ 4 9i(—PdCos9j 4 p s sin 6j)} (4.7) 

(4.8) 

4.3 VELOCITY AT POINT ’P’ 

The absolute velocity vector V, at point ’p’ on the deformed beam can be written as: 

V = Vh 4 r p 4 {to x f* p ) (4.9) 

where Vjj is the rigid body perturbational translation of the hub center Oh which is 
given as: 

Vh = Ql[Rxe x R + RyCyR + Rz^zr] (4-10) 

Transforming all the unit vectors to 4K-system 

Vh — (VH)x^x4k 4 (VH)yiy4k 4 (Vff)*e 2 4* (4-11) 

where 

{V H ) X = QlR x {cos'ii> k +simp k (p s cos9i + Pisin9i)} + 

CllR y {cos'tp k (-P s cos9i - PiSin9j) 4 sin -0*} 4 
QIR Z {P P 4 Pd cos 6 1 - P s sin 9 } } 

( V H ) y = SllR x {cos‘ip k (P s cosQ I + p d sin9i - PpSinBj) - sin ip k cos 9 j} + 


(4.12) 
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tllRy {cos ip k cos 6j + sin ip k ({3 S cos 0/ + sin 6j - (3 P sin 6j } + 

£ }lR z {sin9j} 

(Vh) z = BllR x {cos cos 8; -f fi s sin 6; — /3 P cos Qj) + sin ip k sin 9j) + 

£11 Ry {cos ijjk{— sin#/) + sinV’*(— Aicos0/ + /3 S sm6j — /3 P cos 6 /} + 
£2ZJ? z {cos0/} 

(4 


= e l4 fc/Q[— aOj (—/3 S sin 0/ + (3 d cos 0/ )+cosA ii cosA Q u*] 

+Vk sin A s cosA a — w k sin A a + Tj{—v' k cosA s cos A a cos(0 G + 4>k) 

+v' k 4> k cos A s cosA a sin {6g + <Pk) ~ w' k cos A s cosA a sin(0 G + 4> k ) 

—w' k cosA 5 cos A a <j> k cos(0 G + Ok) + 4>k sinA s sin(0 G + (fi k ) 

+(j> k cos A s sinA a cos (9 a + <p k ) } 

+C {—w'k cosA s cosA a cos(0 G + <^) +w k cosA s cosA a sin (9 G + 4> k ) 
-\-v' k cos A s cosA a sin(0 G + <t>k) + v' k cosA s cosA a (j> k cos(0 G + 4> k ) 

+ sin A s <p k cos(0 G + 4> k ) ~4>k cos A s sinA a sin (0 G + <£*) }] 

+e y4 kl£l[-a6i sin Qj + eidj (—/3 S sin 0/ + /3 d cos 0/ — fi p cos 0/ ) 

-t-e 2 0/ {—fis sin 9 1 + Pd cos 6j )-u k sin A s 4- v k cos A s 

+r){-v' k sin A s cosA Q cos(0 G 4- <£*) + v' k sin A s cosA a <f> k sm(0 G + <p k ) 

-w'k sinA s cosA a sin (9 G + <f> k ) - w' k <j> k sinA s cosA a cos(6 g 4- <f> k ) 

-(j) k cosA s sin (0 G + <f>k) + 4>k sinA s sinA a cos(9 g 4- <f> k ) } 

+({-u4. sinA s cosA a cos {9 G + 6 k ) +w' k j> k sinA s cosA a sin(0 G 4- 6 k ) 
+v' k sin A s cos A a sin {9 G + <f>k) 4-t’* sinA s cosA a <j>k cos(0 G 4- 4> k ) 

-fa cos A, sin (0 G + <f>k) ~ 4>k sin A, sin A a sin(0 G + <p k ) }] 

+e zik lfl[-a9 r cos 6] +ei6j (B d sin 0/ + J3 S cos 6j 4- fi p sin 0/ ) 
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+e 2 #/ (Pd sin &i + P s cos 0/ ) + u k cos A 2 sinA a 
+Vk sin A s sinA a + w k cos A a 

sinA a cos(0 G + <£*) -v k fa sinA a sin(0 G + -hit;* sin A a sinf^c? -H 0*:) 
+w' k 4> k sinAa cos (& G + fa) + fa cosA a cos (0 G + fa) } 

+C{^ sinAa cos (9 G + fa) ~w' k sin A a fa sin (9 G + fa) - v k sinA a sin(0 G + <£*) 
-v' k sin A a fa cos (0 G + fa) -fa cosA a sin(0 G + fa) }] 

(4.15) 


and 


(u x r) = 


&xAk €y4k &zAk 


Ux 
r X 


u v 


r z 


(4.16) 


where r x , r y , r 2 and ui x , u> y , u z are the x,y and z components of f p and u respectively. 
Substituting various quantities in Eq. 4.9 from Eqs. 4.2 ,4. 6-4. 8 and 4.11-4.16, the 
velocity at point p is obtained. This can be written in symbolic form as: 


V = ftl[V x e X 4 k + V y e y 4 k + V 2 e 2 4fc] 


(4-17) 


The components of the velocity vector are given below. 
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v x = {V H )x + (fp)x + (Wj,^ - w z r y ) 

= .R x {cosV>* 4- simp k (/3 s cos6] 4- j3 d sin9i)} 

+Ry{cosipk(— P s cos6i — fidsinOj) + sim/j k } 4- 
R z {Pp 4- PdCosOi - p s sin9i} 

+adj(/3 s sin0i — fidCosOj) + u k cosA s cosA a 4- v k sinA s cosA a — w k sinA a 
+T][cosA s cosA a {—v' k cos(9 G 4- (j> k ) 

+v k <j) k sin(9 G + 4>k ~ w k sin(9 G 4- <j> k ) 

- w' k (j) k cos{0 G 4- 4> k )} + <j> k sin(9 G 4- <f> k )sinA s 
+4> k cos(9 G 4- (j) k )cosA s sinA a ) 

+C[cosA s cosA a {-w' k cos(6 G 4- <j> k ) 

4 -w k <p k sin(6 G + <j> k 4- ^sm(0 G 4- 4>k) 

4 -v' k (t> k cos(6 G 4- <f) k )} + <t>kCOs(9 G 4- <t>k)sinA s 
-<j> k sin(9 G 4- <t> k )cosA s sinA a ] + ( 9j((3 s cos9 j 4- /3 d sin9 /) 

4 -cos9i(—9 x sinipk + 9 y cosip k ) + 9 z sin9j) 

[—asin9j 4- (e\ 4- e<i ){—fidCOs9i 4- /3 s sin0j} 

—eiP p cos9j 4- ( Xk 4- u k )cosA s sinA a 4- v k sinA s sinA a 
+w k cosA a + r)[v' k cos(9 G 4- <j> k )sinA a 
4 -w' k sin(9 G 4- <j) k )sinA a 4- sin(9 G 4- (j> k )cosA a ] 

4-C [w' k cos(9 G 4- (j>k)sinA a - v' k sin(9 G 4- 4> k )sinA a 

+cos(9 g 4- <^*)cosA a ]] — (theta j{-(3dCos9j 

+/3 s sin9i ) 4- cos9j( 1 4- 9 Z ) 4- sin6i(9 x simj> k — 9costJ) k )) 

[acos9j 4- (ej 4- e 2 ){/?iCos0; 4- P d sin9i} 

-eiP p sin9j - (x k + u k )sinA s 4- v k cosA s 
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+7][sinA s cosA a {— v' k cos(6a + (t>k ) ~ w' k sin(6G + <^t)} 
+cos(6c + <f>k)cosA 3 + sin{6c + <j>kStnA s sinA a ] 

+C [sin A s cos A 0 {— tUfc cos [Oq + <f>k ) + v' k sin(OG + (j>k)} 
—sin(6c + 4>k)cosA s + cosific + (f>kSinA s sinA a ]] 


( 4 . 18 ) 
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“I" (^p)y "I" { w x r z fL 'z r x) 

= R x {cosipk (As cos#/ + Ad sin#/ — Ap sin#/ ) — sin 0* cos#/ } 

+jRj,{cos?/>fc cos#/ -fsinti'* (,3 S cos#/ + Ad sin#/ - 6 P sin#/ )} 

+i? z sin#/ -a#/ sin#/ +#/ {(e 1 +e 2 )(-As sin#/ -f Ad cos#/) 

- A P e/ cos #/ } - u fc sin A s + u* cos A, + tj [sin A* cos A a 
{-v' k cos (#c + Afc) + 4 <4 sin(#c 4- Afc) - w'k sin (9 G + <f>k) 

-w' k (j > * cos(#c + A fc ) }-Afc cos A, sin(# G 4- <t>k) + 4>k sinA 4 sinA a cos(# G + <f>k) ] 
+C[sinA, cosA a {-«4 cos (# G + Afc) + w' k 4> k sin(# G + <j> k ) 
v'k sin(# G + (f> k ) + v k d> k cos (# G 4- (j> k ) } 

-<j) k cosA 4 cos(# G + A fc ) - fa sin A, sinA a sin(# G + 0 fc ) ] 

4- (#/ {-Ad cos#/ +A» sin#/ } + (l + # z )cos#/ +sin#/ {# x sin^* - # y cos Ait }) 

71 — 1 

[— a(As cos#/ 4- Ad sin#/ ) + + e 2 + ^ 

e=l 

+(a:jfc +Ufe)cosA s cosA a +y fc sinA 4 cosA a -w fc sinA a 
+77[cosA 4 cosA a {-y* cos(# G 4 <j> k ) — w' k sin (# G 4- <4) } 

- cos(# G + 4> k ) sin A, + sin(# G + <f> k ) cos A s sin A a ] 

+C[cosA 4 cosA a {-*4 cos(# G 4- (f> k ) +v k sin(# G + 4> k )} 

+ sin A 3 sin(# G + <4) +cosA s sinA a cos(# G 4- <f> k ) ]] 

- (#/ - (A P + Ad cos #/ - A* sin #/} — #i cos ipk - # y sin 0* ) 

[—a sin#/ + (ex 4- e 2 ){-Ad cos#/ 4- As sin#/} — ex Ap sin#/ 

+(xjt + m* : )cosA j sinA a 4- v k sin A s sinA a + w k cos A a 

4-77(4. cos(# G 4- (f> k ) sin A a 4 - w' k sin(# G 4- <f > k ) sin A a + sin(# G 4- <p k ) cos A a ] 

+C [~ v k sin(# G 4- <f> k ) sinA a + w k cos(# G 4- 4> k ) sinA a 4- cos(# G + 4) cos A tt ]{4.19) 
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v z = ( V H ) Z + (r p ) z + (w x r y - w y r x ) 

= Rxlcosipk {-p d cos #/ + /3 a sin #/ - P p cos #/ ) -f sin ip k sin #/ } 

+Ry {cos ipk - sin 6j + sin xp k (~/3 d cos #/ -f p s sin 9, - (3 P sin#/ )} 

+4 cos#/ —a#/ cos#/ +#/ {(ei + e 2 )(Ai sin#/ + /? s cos#/ ) 

Pp ei cos #/ } + u k cos A s sin A a + ik sin A s sin A a + Wk cos A a + 7?[sin A 0 
{v'k cos(# G + (t) k ) - v' k (t>k sin(#G + 4>k ) - w' k sin(# G + <j> k ) 

+w[ <t>k cos(#G + 4>k) } + 4>k COS A a COS (#G + 4) ] 

+C[sinA a {w* cos(#g + <£*) - ^ <i>k sin (# G + <pjt) 

-ujt sin(#c + <f>k) + v' k 4> k cos(#g + &) } - <t>k cos A a sin(#G + ] 

+ (#/ + {Pp + Pd cos #/ - & sin #/} + #i cos ip k + 9 y sin xp k ) 

[a cos#/ + (P s cos#/ + Pd sin#/ )(ei + 62) — e\P p sin#/ — (x k + u k ) sin A s 
+n fc cosA 3 +77[sinA s cosA a {-v' k cos (#g + <£*) - u;* sin(#G-f0*) } 

+ cos(#g + <^fc) cosA s + sin(#c + <p k ) sin A s sin A c ] 

+C[sin A s cos A Q {-to* cos(# G + (f>k) +v' k sin (# G + <£*)} 

-sinA s cos(#g + 4>k) + sin A s sinA a cos(# G + 4>k) ]] 

— {9 1 {P s cos#/ 4- p d sin#/ } — (# x sin xp k — 0 y cos xp k )cos #/ + 8 Z sin#/ ) 

n~\ 

[-a{/3.s cos#/ + & sin#/ } + (ei + e 2 ) + 52 ^ 

e=l 

+(xjt + u*)cos A s cosA a +^sinA 4 cosA Q ~w k smA 0 
4- 77 [cos A s cosA a {~v' k cos(# G 4 4>k) ~ w k sin (#g 4 </>*) } 

- cos(#g + 4>k) sin A s 4- cos A s sin A a sin(#c + <f > k ) ] 

+C[cosA s cosA a {+U* sin(#G 4- 4>k) -u-' k cos(#g + <j> k ) } 

4- sin(#c + 4 > k ) sin A s 4- cos A s sin A a cos(# G 4- <£*) ]] 


(4.20) 



Chapter 5 

EQUATIONS OF MOTION FOR 
ROTOR BLADE 


The coupled equations of motion of the rotor blade can be derived using Hamilton’s 
principle. The mathematical form of Hamilton’s principle is stated as follows: 

f 2 {6U - 6T - 6W e )dt = 0 (5.1) 

Jtl 

where U is the strain energy; T is kinetic energy; W e is the work done by the non- 
conservative external loads. 

In this chapter, the expressions for the variation of kinetic energy and strain energy 
of the rotor blade are derived. 


5.1 KINETIC ENERGY OF BLADE 

The total kinetic energy of the beam, T, is defined as : 

r = j Zjfl.ov-vdvoi 

= 5?jT7 S A ^- 9iriiiix 


(5.2) 
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where V is the velocity of an arbitrary point ’p’ on the blade cross section with 
respect to the inertial reference frame. The variation of kinetic energy', 5T can be 
written as: 

ST = VCfL pV.SVdrjdCdx (5.3) 

— ♦ 

Substituting for the velocity V from Eq. 4.17 and integrating ST by parts with respect 
to time, yields 


ST =zL 




J J p[Z u 5u + Z v 5v + Z w Sw + Z' v Sv' + Z’JSw' 4- Z^54>}drj dx (5.4) 


where the terms Z u , Z v , Z w , Z' v , Z' w , Z <p are the coefficients of 5u, Sv, Sw, Sv 1 , 5w', Sep in 
the variation of kinetic energy expression. 

Integration of the expression over the cross-section yields: 


ST = mQ. 2 l 3 ?f [Z u Su H- Z v Sv T Z w Sw + Z' v Sv' -r Z'^Sw 1 4- Z^Scp] dx (5.5) 

After eliminating the higher order terms, the expressions for Z u ,Z v ,Z w ,Z' v ,Z' w ,Z <i> 
are obtained. 
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5.2 STRAIN ENERGY OF BLADE 

The formulation of strain energy outlined in this section essentially follows the prece- 
dure given in Refs. 24 and 28 .// 

5.2.1 Strain Energy 


The strain energy of a beam element is given by 




Crx 

T 

&XX 

* 7xi 7 | 

> i 

| °XT] ' 

. 7xc J 


l . 


> dr) d( dx 


5.2.2 Explicit Strain-Displacement Relations 


(5.6) 


The expressions for non zero strain components written in terms of u,v,w and <p are 
given as [24, 27] 

= U X + + ^(?? 2 + C )4>1 + 

+ar 0 ( ~ - [vcos (6 g + <f>) - £ sm(6 g + 4>)]v xx 

-[? 7 sin( 0 s + 4>) + C cos(0 fl 4- <f>)]w xx + 77 ( 7 ^, x - t 0 %^) + C(7xc,x + r 0 j xr) ) 


7xi 7 — 7xtj "h Oc^ff C(^*x “h 4*o) 
7xC = 7xc + - V{<t>x + 4>o) 


where 


4> 0 = {v xx cos 6 g + il'xt sin 9g)(-v x sin 6 C + w x cos 6 G ) 


The underlined term in e xx represents the axial strain at the elastic axis. These strain 
expressions can be simplified using the following assumptions: 
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• The transverse shear at the elastic axis is assumed to be zero. 


The wraping amplitude a is assumed to be equal to -<f> x 


The simplified strain components can be written as 


e*x = u x + )jvl + ^ w 2 x + ^(v 2 + C 2 )4>l - ^<f>xx - - ^c)]^ 

[77 cos(0 G + (j) k ) - c sin (0 G + <pk) ]v xx ~ [r] sin(0 G + <j> k ) + C cos(0 G + <f> k ) ]w xx 

yxr] — "h O^x Q(f>o 

7x< = -(^c + V)<Px + V<l ) o 

<t >0 = ( v xx cos 0 G + w xx sin 6g){-v x sin 0 G + w x cos 0 G ) 


5.2.3 Stress-Strain Relations 


Assuming that the blade is made of isotropic material. The stress-strain relationship 
is given by the following equations: 

1 \E 0 0 H e M 1 

(5.7) 


O xx 


' E 0 O' 


Zxx 1 

< &XT) | 

► — 

0 G 0 

0 0 G 

< 

I 7X7? j 

l 7xc J 


5.2.4 Strain Energy Variation 


The variation of strain energy of the i th element is given by 

r '\ T 

w -W// 


$€xx 

. hxc ) 

The variation of the strain components is given as: 


'-'XT] 

OxC, 


dr] d(, dx 


(5.8) 



28 


5e xx = 6u x + v x 6v x + w x 5w x + (rf + C 2 )0x^n — ^5<j> xx - [r 0 (C^^ — rjj^/ ^)]5<j> x 
- (rj cos 0 G - Csin 0 G )(Sv xx + <t>5w xx + w xx 8<t>) 

— (rj sin 9 G + Ceos 9 G ){8w xx - </>8v xx - v xx 8b) 

hxT) = -(^ + C)^X“C^O 

<HcC = ~{^c~v)6(t> x + r)5<j> 0 

where 


8(p 0 = {-8v x sin 0 G + 6w x cos d G ) {v xx cos 9 G + w xx sin d G ) 
+(— u x sin 6 g + w x cos 9 G )(8v xx cos 9 g + 8w xx sin 9 G ) 


It is assumed that the variations of strain components are of the same order as the 
corresponding strain components. Substituting the above expressions in the strain 
energy variation and integrating over the cross-section, the variation in strain energy 
for a beam element is obtained. The expression for 8b\, interms of stress and moment 
resultants, is given as: 

r lc - 

8U t = E 0 l I {V x (8u x + v x 8v x + u: x 8w x ) 

Jo 

+(—t 0 P x - R x + S x + T x <p x )8(p x + S'Mo ~ Px^bxx 
+[v xx {M!j cos 9 g - Ml- sin 9 G ) -f w xx (M' v sin 9 C + cos 9 G )]5<f> 
+[{Mt] sin 9 g -I- Mq cos 9 G ) + 0(Af' cos 9 G — sin 9 G )]8v xx 

4-[(-Mt; cos 9 g + Me, sin 9 G ) + <£(M' sin 9 G + cos 9 G )]8w xx }dx 



Stress and Moment Resultants are defined as: 


V x = J J {Ee xx )dr,dQ 
= EA[u x + ^v 2 x + ^wl] 

+±EAC 0 [<£) - EAD 0 [<j) xx ] - ToEAD'M 
-(EArf a - <j>EAQ\v xx ] - (EAQ + <f>EAr] a ){w xx } 

(S 

Mr, = J J (CEe xx )drjd<: 

= EACa[u x + ^v 2 x + ^wl] 

+^EAC 2 [<t>l]-EAD 2 [(f> xx }-T 0 EAD' 2 [(f> x ) 

- (EI vC - <pEI m ) [v xx ] - (EJ m + 4>EI«) W 

( 5.1 

m; = //«•*„« 

= EAC a [u x + ^vl + ^w 2 x ] 

- EAD 2 [<f > xx ] - T 0 EAD' 2 [(t> x ] 

{EJfjQ ) ['Un] (EI r j J ]')[w xx \ 

( 5 . 12 ) 

M c = J J ( ~TjEe xx )dr]d ( 

= -EAv a [u x + h'l + ~w 2 x } 

-\eACMI) + EAD 1 [d> xx ] + ToEA&M 
+(£/« — + (El( 1 , + Ill'll] 


( 5 - 13 ) 
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= / / (-T)Ee xx )drid(; 

+^A[^x] + r 0 ^D;[^] 

+ (^Cc)bn] + (-E^CnH^xx] 

P x = J J (<HEe xx )dr]dC 

= +EAD 0 [u x + l -vl + \wl ] 

+ijBAD 4 [^] - EAD,[&J - r 0 EAD b [(j> x ) 

-{EADl - <j>EAD 2 )[v xx ) - (EAD 2 + <j>EAD[)[w xx } 

K = J /[(C*, ~ fi* ( )Ec„))dndC 

= +EAD' 0 [u x + l -v\ + ^wl] 

+ l -EAD\[<j>l \ - EAD b [<t> xx ] - ToEAD'M 

-(. EADi - <£EAD~ 2 )[v xx ) - (EADl + <fEAD~()[w xx } 

T x = J j [(rf - C 2 )Ee xx )]drjd( 

= +E AC 0 K + ^ v 2 x + ^ w 2 ] 

a\eAC z [ 4>1] - EAD 4 \<j> xx ] - r 0 EAD'S x ] 

~{EAC{ - <j>EACl)\v xx \ - (EACl + <fEAC()\w xx ) 

R x = J J - *c7 *c)G]<MC 

= (Gh = GJ 2 )[4> x \ + GJM 


(5.14) 


(5.15) 


(5.16) 


(517) 


(5.18) 
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Sx = J J[{vixc ~ C7«j)G]^dC 

= (GJo = GJ, )[<!>:] + GJM 

( 5 . 19 ) 

S'x = j J [(VJxc - Clxr,)G}dridC 

= (GJ 0 = GJJfa) 

( 5 . 20 ) 

In the above expressions, the moment resultants M', Ml and S' x have the same defini- 
tions as M v ,Mq and S x respectively. But higher order terms are neglected in defining 
M^ . and S' x , since these expressions are coupled with terms of order e 2 . Whereas 
Mr),M^ and S x are coupled with terms of order e. The cross-sectional integrals asso- 
ciated with the strain energy variation are defined as: 

EA = j f[E}drjd<: 

EAr] a = j j[E V }dr 1 d C 
EAQ = / / [EQdrjdC 
El m = j flE^dvdC 

EI« = J fiEtfidridC 

EI v q = / / [EriQdrjdC 
EAC 0 = f J[E( V 7 + e)\drid( 

EACi = j f[Er 1 (ri 2 + C 2 )}dvd(; 

EAC 2 = J f[Eav 2 + ?)}drid<; 

EAC 3 = J J[E(t ] 2 + c 2 ) 2 3^C 
EAD 0 = J JlEV}d V d<; 

EADi = J J[Eij&]dr]dQ 



ead 2 

ead 3 

ead 4 

ead 5 

EAD ' 0 

EAD\ 

EAD ' 2 

EAD' Z 

EAD' a 

GJ 0 

GJi 

GJ 2 


= J jlEWdndt, 

= J J [E* 2 ]dr)d( 

- J J i EI ' r i'‘ + ?)md< 

= J /[£(<«'„ - 

= J fm% - ^ ( )}dr,i( 

= / /[£C(C*, - »4 C )]^C 
= / /[£(C*, - 

= J J{E(V 2 + ( 2 )(( < S'n- r )V!«a)}drld( 
= J /lG(i ) 2 + ( 2 )]*)<iC 
= / /[G((», - i)*<)]*XiC 
= j J[G^l + ^\dr,ii 


Certain additional terms are defined as: 


£j4r? a = EAri a cos $g — EM£ a sin d G 
EAQ a = £-4% sin 8 G 4- EAC, a cos 0 G 
EACl = EACi cos $g - £MC 2 sin $g 
EAQi = JEMCi sin $g + £-402 cos #g 
EADi = £MZ?i cos 0 g - EAD 2 sin 8 G 
WAD 2 = EAD\ sin #g + EAD 2 cos #g 
EAD[ = EAD\ cos 0 G - si n 
WAD l = EAD\ sin 8 q + EAD 2 cos 6 q 
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El* 

= El* 

cos 8q 

- Elm sin 


Elm 

= Eire 

sin Oq 

"f* Elfjfj cos 

@G 

El* 

= EI« 

cos 9q 

— El ^ sin 

e G 

El c. 

= £/« 

sin 6q 

+ EI{ V cos 

@G 

GJ 

= GJ 0 

-2GJ 

i + gj 2 





Chapter 6 


FORMULATION OF ELEMENT 
MATRICES ASSOCIATED WITH 
KINETIC AND STRAIN 
ENERGY VARIATION 

6.1 FINITE ELEMENT DISCRETIZATION 

The variational expressions associated with the kinetic and potential energy of the 
rotor blades are nonlinear. The unknowns are the deformation functions u,v,w, and <j>. 
These are dependent on both space and time. The spatial dependence is eliminated 
using a Rayleigh Ritz finite element formulation. The blade is divided into subregions 
(finite elements) as shown in Fig. 10, and the total dynamic potential is calculated 
for each subregion. By applying Hamilton’s principle to each subregion, a discretized 
form of the equations of motion can be obtained. In this development, each subregion 
is modelled by a straight beam type finite element. These beam elements are located 
along the elastic axis (line of shear centers) of the blade. 
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The discretized form of Hamilton’s principle is written as: 

ft*" 

/ Y,{6Ui - 6Ti - SW ei )dt = 0 (6.1) 

Jt i »=i 

where: 

N = Total number of finite elements in the model. 

5Ui = Variation of the strain energy in the i th element. 

6Ti = Variation of the kinetic energy in the i th element. 

6W ei — Virtual work of the external loads in i th element. 

Assume that the four unknown displacement functions of the beam finite element are 
expressed in the following form. 


V 


■ {<M T o o o ' 


{V'} 1 

w 


0 {K} T 0 0 


on 



o 

■e- 

v -^- f 

o 

o 


w 

u 


1 

t-1 

a 

O 

o 

o 

1 


L m J 


where {<£„}, {<p w }, {<^,}, {<£u} are space dependent interpolation functions; and 
{V}, {IV}, {$}, {U} are time dependent nodal values of v, w, 6. u, respectively. 


(6.2) 



r Vi 

l 


'. W\ 


f <h I 


f «1 1 

m = < 

V 1 

V2 

. v' 2 , 

. ; {W} = < 

W 1 

W2 

, n . 

«;{*} = 

l * J 

■;{!/) = < 

u 2 > 
l U 3 , 


The nodal coordinates are shown in Fig.ll. 

The variation of the displacement functions for the beam can be written as: 


Sv 


' {<fi v } T 0 0 0 ' 


<«'} 1 

6w 


o {<?U T o o 


{««•} 

6(1) 


o o {^} T o 


{»} 

6u 


1 

bl 

r-H 

o 

o 

o 


{« ! } J 


(6.3) 


In this development, Hermite interpolation polynomials are used for the space de- 
pendent interpolation functions. A cubic Hermite interpolation polynomial, {$ c }, is 
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used for the bending deflections (v,w) and a quadratic Hermite interpolation polyno- 
mial, {$ g }, is used for the torsional rotation and the axial deflection (u). The 


mathematical expressions for these polynomials are given as: 


v } — 


1 - 3£ 2 4 - 2£ 3 

ut-2e+e) 

3£ 2 - 2 f 3 

u - e + e ) 




where: 


{4><t>) = {<t>u} = 


- 3£ + 2£ 2 

ii-Ae 

- z+e 




(6.4) 


(6.5) 


£ — %k/le 

Xk = spanwise (axial) coordinate of the beam element . 
l e = length of beam element. 

For bending deformations, the nodal parameters are the displacements and slopes at 
both ends of the beam element. Therefore, the resulting elements will have interele- 
ment continuity for both displacements and slopes. In addition, because of the cubic 
interpolation polynomial, bending strains vary linearly over the element length. 

The quadratic Hermite interpolation polynomial is used for torsional rotation (0) and 
the deflection (u). This polynomial has the capability of modelling a linear variation 
of strains along the element length and therefore provides the same lavel of accuracy 
as the beam bending element. The nodal parameters for these elements are chosen as 
the values of the displacements function at the two end nodes and at the mid-point 
of the element. 

The resulting beam element has 14 degree of freedom: 4 in-plane (lag) bending de- 
grees of freedom, 4 out-plane (flap) degrees of freedom, and 3 degrees of freedom each 
of torsion (<£), and axial deflection (u). The nodal degrees of freedom are shown in 
Fig.ll. 
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6.2 ELEMENT MATRICES ASSOCIATED WITH 
KINETIC ENERGY VARIATION 

The beam element matrices associated with the kinetic energy variation are derived 
by substituting the assumed expressions for the displacements funtions (Eq. 6.2) in 
the kinetic energy variation ST (Eq. 5.5) and carrying out the intergration over the 
length of the beam element. The resulting variation of the kinetic energy can be 
written in the form: 


ST, = -{<Ss} t ([m]h*h{?} + [m c ]„ x11 «} + [a:'T x11 {9} + {V 1 },,*! 


[R x } 


R x } 


f ?x] 


e x 

i 

► + [M 2 ]ux 3 < 

*• 

{ R z J 

► + [A/ 3 ]i4x3 * 

fj 

> + [A^ 4 ]i4x3 < 

9y 

{ 9z J 


+ 0 7 }l4xl + {V’ JVi }i4xl) 

(6.6) 

Where {q} represents the vector of unknown nodal degrees of freedom 


I?} 14x1 



(6.7) 


Detailed expressions for the various matrices defined in Eq. 6.6 are given as follows: 
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6.2.1 Mass Matrix [M] 14xl4 

[Mu] = J [(cos 2 A s + sin A a cosA a + sinA 5 sin A a )m{<3> c } {$ C } T 

— (sin A„ cos A a cos A s + cos A s sin A s cos 2 A a — sin A 0 sin A s cos A a ) 
{mr) m cos(e G + (f> k ) - mC m srn(e G + </>*)} {$ c } {$c} T 
— (cosA s cos 2 A a sinA s +sinA s cosA 0 cosA s — sinA s sin 2 A a ) 

{mr} m cos(0 G + 4> k ) - m sin(0 G + <f> k )} {$[.} {$ C } T 
+(cos 2 A s cos 2 A a + sin 2 A 5 cos 2 A a +sin 2 A Q ) 

({/m cc cos 2 (0 G + <j> k ) + Im^ sin 2 (0 G 4- <p k )} 

-{Inn ~ hi) sin(0 G + <j>k) cos (0 G + <f> k ) ){$' c } {$'} T ]dx 

[M 12 ] = f [(— cos A* sinA 5 cosA a — sin A* cos 2 A a cos A* +sinA 5 sin 2 A a ) 

Jo 

{ra77 m sin(0 G + <t> k ) + rnCm cos(0 G + <f > k )} {$ c } {$'} r 
+(cosA s cosA a sinA a + sinA a cosA a ) 

{mq m cos(9 g + 4>k) - mCm sin(0 G + <£*)} {$' c } {$ C } T 
+(cos 2 A s cos 2 A a +sin 2 A s cos 2 A a +sin 2 A fl ) 

{{I m - /cc) sin(0 G + <f>k) cos{6 g + <j> k ) + 7„ c {cos 2 (0 G + <j> k ) - sin 2 (0 G + 4> k )} ) 

rU 

[M13] = / [(- cos 2 A s + sin 2 A s cos A a ) 

Jo 

{m7? m sin(0 G + (j> k ) + m£ m cos(0 G + 4> k )} {$ c } {$ 9 } r 

+(cosA s sinA s sinA a 4- sin A s cosA a cosA s sinA a +sinA s sinA a cosA 0 ) 

{m77 m cos(0 G + <f > *) — m(mSin(0 G + $*)} {$ c } {$ 9 } r 

(- cos 2 A s sin A a cos A a - sin 2 A s cos A a sin A a 4- cos A a sin A a ) 

({/m^ cos 2 (0 G + 4> k ) 4 - lrn m sin 2 (0 G + <f> k )} + {2 I v q cos [0 G + <f> k ) sin(0 G + &)} ) 



39 


TO to t ]dx 

rU 

[Mu] = / [(— cosA s sin A* + sinA* cosA* cos 2 A a 

Jo 

+ cosA s sin 2 A a sinA s )m{$ c } {$ 9 } r 

+(— cos 2 A 5 cos 2 A a +sin 2 A s cosA a + sinA* sinA a cosA a ) 

{mrim cos {0 G + fa) - mCm sin (0 G + fa)} {$[.} {$J r }dx 

[M 21 ] = [Mu} T 

[A/ 22 ] = [(cos 2 A a - sin A s sin A a )m{$ c } {$ c } r 
(cos A a sin A a + sin A s cos A s cos A a ) 

{mrjmCos(O c + fa) - mC m sin{0 G + fa)} {$ c } {$c} r 
{mrim sm(0 G + fa) + mQ m cos (0 G + fa)} {$ c } {$' C } T 
(cos 2 A s cos 2 A a 4 - sin 2 A s cos 2 A a + sin 2 A a ) 

{{Im^ sin 2 {6 g + fa) + Im vr] cos 2 (6 C + fa)} + { 21 ^ cos{0 G + fa) sin (6 G + <£*)} ) 

m TO 7, )dx 

rh 

[A/ 23 ] = / (cos 2 A a — sin A s sin A a cos A s ) 

Jo 

{mT ] m cos (0 G + fa) - mQ m sin( 0 G + fa)} {$c} {$ 9 } X 
- sin 2 A s {mri m sm(0 G + fa) + m cos(0 G + fa)} {$ c } PQt 
(cos 2 A s cos A a sin A 0 + sin 2 A s sin A a cos A a - sin A a cos A a ) 

({I m - Itf} sin(0 G + <Pk) cos(0 G + fa) + 1 v ^{cos 2 {0 g + <j> k ) - sm 2 (0 G + <?*)} ) 

{$c} {^?} r } dx 

ft/ 24 ] = f‘ [(cos A* cosA a sinA a -sin A, cos A* cosA 0 ) 

Jo 



m{$ c } {$ 9 } r 4-(-cos 2 A 4 cosA a +sin 2 A 4 cosA a 4-sin 2 A a cosA 4 ) 
{mrjm sin(d G 4- <f>k) + mCm cos (0 G + fa )} {$' } {$J T ]dx 

[M 31 ] = [M n ] T 

[M 32 ] = [M 23 ] 7, 

/*Z C 

[M 33 ] = / [(cos 2 A a 4- cos 2 A s sin 2 A a — sin 2 A s sin 2 A a ) 

J 0 

"+■ ^vv)i^g} {^g } 7 

( 2 sin A s sin A a cosA s )(/ 7 , c {cos 2 (0 G + <t>k) ~ sin 2 (0 G 4- fa)} 

~{Ir)T) ~ -^cc} Sin(^G + fa) COS(0 G + (j) k) H^g} {^g} 7 ^ ]^ T 

[A/ 34 ] = [ [(cos 2 A 4 cosA a sinA 4 4 -sin 2 A 5 cosA a cosA 5 ) 

Jo 

{mT) m sin(0 G + <j>k) + mu cos(0 G 4- fa )} {$J {$ 9 } T 
4 - (cos A s cos A a sin 2 A a 4- cos 2 A s cosA a sinA a ) 

{mrj m cos (6 g 4- fa ) - mCm sin(0 G 4- 0*)} {^ 9 } T ]dx 

1 M 41 ] = [M 14 f 
|M 42 ] = [M 24 ] t 
|M 43 ] = [M^f 

(M443 = f l ‘ (cos 2 A cos 2 A a - sin 2 A s 4- cos 2 A, sin A, sinA 0 )m{$ g } {$J r 
Jo 
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6.2.2 Matrix [M c ]u x u 

fit 

[iV/fi] = / m[— sin A 3 cosA a cos A, -cos A, sin A, cosA a 

J o 

+ sin A, sinA a cos A, ]{$ c } {<3> c } r 

+[—mr) m cos(0 G + <f> k ) + mCm sin {9 G + <j> k ) ]{cosA, cosA a 
+ sinA, sinA a cosA a }{$ c } {$' C } T dx 

fit 

[A/jy = / m cosA s cosA a {$ c } {$ c } r dx 

J o 

[ jV/ l 3 l = [ [2 cos A, cosAa {Irm +/ C c ) sin (^G + 0Jfc) cos(^ G + 0fc) { $ c} {^?} T 

J 0 

-cosA a cos A, Itf (cos 2 ((fc + &) - sin 2 (ffa + $k) )]dx 

[A/fJ = 0 

fU 

[A/ 2 C i3 = / [cos A s cosAa {(rnrj m cos(9 G + <p k ) ~ sin [9 G + (j> k ) ) 

J 0 

-{mrj m sin (9 G + <p k ) + m cos(0 G + <f>k) )}((A cos 9j - /? d sin 9i + 1) ){$ c } {$'} r 
4- cosA a cos Aj {mr) m sin {0 G +<l>k) -mQ m cos (9 G + 4> k ) }{<£<;} {$' C } T ]dx 

rl t 

[M% 2 ] = / [cos A, cosA 0 {mTj m sm(9 G + (j> k ) +m(^ m cos {9 G + (f> k ) } 

Jo 

((A cos 8i - /3d sin 9i + 1) ) {$ c } {$ C } T ]dx 

[M| 3 ] - [‘ [-{mrjm sin(0 G + <f> k ) + mCm cos (9 G + <j> k ) } cosA a {$ c } {$q} T 

J 0 
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-sin A, sinA a {mr} m sin(0 G + (j> k ) +mQ m cos {9 G + fa) }{$' c } {$J T 

— cos A, sin 2 A a {mrj m sin (d G + <f> k ) 4- rn( m cos(0 G + fa) }{$' } {$ q } T 

— cos A, cos A a (/^ — Ijp j )(sin 2 (0 G + <f> k ) — cos 2 (9 G + <f > k ) ){$c} 

+{(^CC ~ Un ~ 21* }sin(0 G + <£*) cos(0 G 4- <£*) cos 2 A, cos 2 A a {$(.} {$g} r ) dx 

[M 2 c 4 ] = 0 

rU 

[ M n] = / cos A, cosA a [{/ <c cos(0 G 4- <£ fc ) - J* sin(9 G + fa) } 

J 0 

(-sin {0 G + fa) cos A s 4- cos (0 G + $ k ) sin A, sin A 0 ) 4- {-/^ cos(0 G + <fo) 

+/ w sin(0 G + ^fc) }(-cos A s cos(0 G 4- <£*) - sin(0 G 4- <£*) sinA 3 sinAa )] 

— cos 2 A 4 cosA 0 [{/ cc - I m }cos {9 G + fa) sin(0 G + fa) 

+/* (cos(0 G + fa ) 2 - sin 2 (0 G + <£*) )]{$,} {$' c } r dx 

[M 3 c 2 ] = [ - sinAa [{/« sin(0 G + 0 fc ) 

4-/* cos(0 G 4- <£*) }(sin(0 G 4- <£*) sin A a +cosA a sinA a cos(0 G 4- 0*) ) 

+{/* sin(0 G + 0*) + /jjjj cos(0 G + fa) } 

(cos(0 G + fa) sin A s - cos A s sin A a sin(0 G 4- fa) )] - ( cos 2 A s cos A a 
4- sin 2 A s cosA 0 ){I^ sin 2 (9 G + fa) + I m cos 2 (0 G 4-0*) 

+2 1^ sin (9 G + fa) cos(d G 4- fa) }- cos A, cosA a [{/<;<■ sin (0 G 4- <£*) 

+1* cos {9 G + fa) }(-cos A s sin(0 G 4- 0*) 4-sinA s sinA a cos(0 G + 0 fc ) ) 

+{I m cos(d G + fa) + 1 sin (9 G + fa) }(cos A, cos(0 G + fa) 

-sinA 3 sinAa sin(0 G + fa) )]{$ c } {$c} T dx 

[M 33 ] = [ cos A, [{/^ sin(# G + 4> k ) + /* cos(d G + fa)} 



43 


(sin(0 G + fa) sinA 4 + cos(0 G + <£*) cos A, sin A a ) -f {/,,£ sin(0G + <^t) 
+/ w cos (0 G + <j>k) }(cos (6 c + <t>k) sinA s -sin(^ G + ^ fe ) cosA 4 sinA a )] 
{*«} {$>J T dx 

f U 

[M 3 c 4 ] = / [-m cos A , cos A a + {mT] m cos (0G + fa) ~ m( m sin {0 G + </>k) } 

J 0 

(cos A s - sin A s sin A a ) - {mr) m sin(0 G + fa) - mC, m cos(0 G + fa) } 
cosA s cos A a ]{$,} {$ g } T dx 

rt e 

[M 41 ] = / m[cos 2 A 4 cosA a + sin 2 A 5 cosA a 

Jo 

— cos 2 A 4 sinA a {^c} T dx 

[M 42 ] = [ m[cos 2 A s cos A a — sin A s sin A a 

Jo 

-cos A a sin A* ]{$«,} {$ C } T dx 

m 3 ] 

m 4 ] 


f - cos A s cosAa ( mr} m cos {6 c + fa) - m( m sin (9 C + fa) ){$?} {$c} 7 dx 
Jo 


rle 9 

/ m[— cos A s sin A a cos A a 

JO 

+ sinA s sinA a cos A s ]{3> g } {$ g } T dx 
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6.2.3 Centrifugal Stiffness Matrix [K $ xl J 

\K n } = f [-m(cos 2 A a cos 2 0/ 

Jo 

4 - cos 2 0 j cos 2 A s sin 2 A 0 

- cos 9] sin 2 A s sin A a cos A a (3 P + /3 d cos 9j - p s sin 9j ) 

{$ c } pet 4- (cos 2 6 ] sin 2 A a cos 2 A a 

4- cos 2 $j cos 2 A s cos 2 A a )({/m a cos 2 (0 G 4 - <?*) 4 - Im m sin 2 (0 G 4 - 0 *)} 
-{2/„c cos(0 G 4- (j> k ) sin (0 G + fa)} ){$' c } {<^} T }dx 

fie 

[K \ 2 ] = / [m(cos 2 6] sinA s sinA a cosA a —cos 6] sinA s cos 2 A Q 

P p 4- PdCOs9j — P s smdj ){$ c } pet 4- (cos 2 9j sinA a cosA A cosA a 
{mr] m cos {9 g + fa) ~ m C m sin(0 G 4- fa)} {$' c } {$' c } r 
-(cos 2 6j cos 2 A s cos 2 A a ){{I VV - 1^} sm(9 G + fa) cos (6 G + fa) 
-/„ c {cos 2 (0 G 4- fa) - sin 2 (0 G + 0*)} ){$' c } {$<} r ]dx 

[iC 13 ] = 0 

[K 14 ] = f [m(cos 2 0 / sinA a cosA 5 sin 2 A a ){$ c } {^ g } T 
Jo 

4 - (cos 2 61 sin 2 A s cosA 0 —cos 2 9] cos 2 A s cos 2 A a ) 

{m 77 m cos( 0 G + fa) - mCm sin( 0 G + fa)} {$c} {$<?} r ]dx 

[K21] = [Ki2] T 

[K22] = [' [-m(cos 2 0 j sin 2 A a sin 2 A s 4 - cos 2 #/ sinA a sinA s ){$ c } {$c} T 

Jo 

4 - (cos 2 9i sin 2 A s cos 2 A a 4- cos 2 9j 
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cos 2 A a cos 2 A 3 )({/m fC si n 2 ( 0 G 4 <t> k ) 4 Im m cos 2 (0 G 4 - fa)} 

+{2/t?c cos(0 g 4- <f>k) sin(0 G 4 (f > k )} ){$(.} {$c} T ] dx 

[K 23 ] = 0 

rl e 

[K 2 4 ] = / [m(cos 2 0 / cosAj cosA a sinA a —cos 0 / cosA 3 

Jo 

cos 2 A a /? p + /? d cos 0/ - sin 0/ ){$ c } {$ 9 } T + (cos 2 0/ sin 2 A, cos A a 

— cos 2 0/ cos 2 A s cos 2 A a — cos 0/ cos 2 A s sin A a cos A a fa + fa cos 0/ - fi s sin 0/ ) 
{mrjm sin(0 G + fa) + m^ m cos(0 G 4- 0*)} {$(.} {$J T ]dx 

[A' 31 ] = [ [{mT] m sin(0 G 4- <£ fc ) 4- m( m cos(0 G 4- fa ) } (cos 2 0/ cos 2 A 4 4- a cos 0/ 

JO 

4-(ei 4- e 2 ){/?i cos 0/ 4- /3<* sin 9 { } - ei/? p sin0/ 

4 - cos 0 / cos A, cosA a {— a(/3 s cos0; 4 -/?dsin 0 / ) 

n— 1 

4 -ei 4- e 2 4- 53 ^ei + £*:}){$?} {$c} r 
1 = 1 

— [(2cos 2 0/ sinA 4 cosA s cos A a )({/„ — /^} sin(0 G 4- fa)cos(9c 4- <£*) 

4-/^{cos 2 (0 G + ^ fc ) - sin 2 (0 G 4- 4)} ) - ({Im cc cos 2 (0 G 4- &) + /m m sin 2 (0 G 4- <£*)} 
+{2/„ c cos(0 G 4 - <£*) sin(0 G + fa)} )(cos 2 0/ cos 2 A, sin A a cos A 0 

- cos 2 0/ cos 2 A a cos A 3 fa 4- fid. cos 0/ - fi s sin 0/ )]{$?} {$c} r ]dx 

[K 32 ] = [ * [(cos 2 0/ sin A, sin A 0 {m? 7 m sin(0 G 4- fa) + mC, m cos (0 G 4- <£k)} 

Jo 

4 -cos 2 0/ sin A s cos A 3 sinA a {mr] m cos(0 G 4- <fo) - mC m sin(0 G 4- fa)} ){$?} {$c} T 
— [(cos0j sin A;, cos 2 A, cos A a )({/,pj — fac} sin(0 G 4- <fi k ) cos(0 G 4- 0*) 

+/ 7 )C {cos 2 (0 G 4 - fa) - sin 2 (0 G + fa)} ) 4 - ({/m« cos 2 (0 G + fa) 4- /m w sin 2 (0 G 4- fa)} 
+{21^ cos(0 G 4- fa) sin(0 G 4 - fa)} )(cos 2 0j cos A, sin A, cosA a 
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- cos 2 9 1 cos 2 A a cos A, fa + fa cos 0, - fa sin 0 r )] {fa} {$^} r ]dx 

rU 

[K 33 ] = y o [(ciit cos 2 A a - cos 2 0/ sin 2 A, sin 2 A a ) 

({Im K cos 2 {9 a + fa) + Im m sin 2 (0 G + 0 fc )} 

-{/m c<: sin 2 (0 G + <p k ) 4- cos 2 (0 G 4 <£ fc )} ){$ ? } {$,} T ]dx 

rU 

[A 34 ] = / [(cos 9 1 sin A 3 cos A, 4 cos 2 0/ cos A s cos A a sin A, ) 

J 0 

{mrjm sin (0 G + <p k ) 4- m(m cos (0 G + <f> k )} 4- (cos 2 0/ cos 2 A s cos A a sin A a ) 

{mr) m cos (0 G 4 fa) - mC m sin(0 G 4 fa)} ]{$,} {$J r dx 

rl e 

[K n ] = / [m(cos 2 0/ sinA 3 cosA s sin 2 A, ){$,} {$ C } T 

J 0 

4 (cos 2 9 1 sin 2 A, cosA a - cos 2 0/ cos 2 A, cos 2 A a ) 

{mr] m cos (9 0 + fa) - mC m sin(0 G 4 fa)} {$,} {$' C } T ]cte 

[A42] = f [m(cos 2 0/ cos A s cosA a sinA a ){$,} {$ c } t +(cos 2 0 j sin 2 A s cosA a 

JO * 

- cos 2 0/ cos 2 A s cos 2 A a ) {mrjm sin (0 G 4 fa) 4- mCm cos(0 G 4- <£*)} {$,} {$' } T ]dx 

[K*] = 0 

[Af 44 ] = [ [-m(cos 2 6 r cos 2 A a +cos 2 0/ sin 2 A, sin 2 A a ){$<,} (<& 9 } r }dx 

Jo 
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6.2.4 Vector [V L ]u x i 

K'n] = f [— { mT] m sin {9 G + (f> k ) + m( m cos(0 G + 4>k) } cos 2 0/ cos 2 \ 3 cos 2 A a 

J o 

{ma(2cos0/ — sin 0/ ) cos A a cos 0/ — m(e i+e 2 )cos 2 A s cos 2 0/ }{$ c } 

— {mrj m sin {9 G + <t> k ) + mC, m cos (0 G + <fo) }cosA, {$ c } 

+{mr/ m cos(0 G + 4>k) - mCn sin (9 G + <fi k ) } cos 2 A, cos A a ( ){$' c } }dx 

[V 2 \] = [ [{-m( ) cos 2 A 3 + (3 d cos 9 1 - p s sin 0/ }{$ c } 

J o 

{ma(2 cos 9[ - sin 0/) sin#/ cosA a {$ c } + {m Tj m sin {9 G + <f> k ) 

+m^ m cos (9 G + <pk) }( ) cos 9 1 cos A, {$'} ]dx 

[V 3 \] = / [{(/« - Ir, v )sin(0 G + <j) k ) cos(0 G +<j) k ) cos 2 9 1 cos 2 A s cos 2 A a 

— I-rtC {sin 2 (#G + 4>k) - COS 2 (#c + 4> k ) } COS 2 0/ COSA a COsA a sin(#G + d>k) 

+I m cos {9 G + (j> k ) } sin 0/ cosA a ]{$ g } dx 

K] = jf [-»*{ + o 

—a( }cos 0/ cos A 3 cos A a {$,} ]dx 



48 


6.2.5 Matrix [M 1 ]i 4 X 3 



f le 

/ m[— cos A s 

Jo 

+ sinA. sinA„ 


cos # j + sin A s cos A a 
sin#/ ]{$c} dx 



K] 

[Aft] 

Aft] 



/•ft 

= / m[cos A s cos 6j - sin A s cos A a 

Jo 

-sinA s sinA a sin#/ ]{$ c } dx 

fle 

= m[cosA s sin#/ +sinA s sinA a cos 6j]{$ c }dx 

fle 

= / m[cosA a sin#/ -sin As ]{$ c } dx 

J o 

fle 

= / ra[-cosA a sin 0; -sin A* ]{$ c } dx 

Jo 

fie 

= / m[cos A a cos - sin A* /3 P + fi d cos 0/ - /?* sin 6j ]{$ c } dx 

Jo 


KM = o 

[Aft] = 0 


[Aft] = 0 

[Aft] ■ jT m [cos A s cos A a +sinA s cos#/ 
+ cosA s sinA s sin#/]{$ 9 } dx 
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[Mi 2 ] = j m[cos A s cos A a - sin A s cos 0/ 

-cosA s sinA a sinflj ]{$ 9 } dx 

[JW«] = f m[cos A s cos A a /? p 4- cos 6j - (3 S sin 6j 4- sin 0/ sin A s 
4-cosA s sinA s cos#/ ]{$ 9 } dx 
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6.2.6 Matrix [M 2 ]i 4x 3 


[m x 2 j] = J m [- 2 cos A a cos 8j - cosA, cos 9j 
+ cos 2 9] sin A, cosA a ]{$ c } dx 


rU 

= / mf-2 cos A* cos 0/ + cos A. cos 8i 

Jo 

- cos 2 8 1 sin A s cosA q - cos A s cos 6j ]{$ c } dx 


f ‘ 

/ m[ cos A s cos 6 1 — sin A s sin A a sin 6j + cos A s cos 6j 
Jo 

+ sinA a cos 6 1 cosA a sin 2 As ]{<E> C } dx 


[Ml] - f: m[— cos A 5 + cos 8 1 cos A a (3 P + /3a cos 8j - 0 S sin 8j 
+ cos 2 dj sinA a + cosA q sin#/ ]{$ c } dx 

[m| 2 ] - f; m[ cos 6j cos A a 4- cosA a 

+ sinA s -sinA 5 ]{$ c } dx 


flc 

/ m[ cos 9 1 sin A a 

Jo 

+ sinA s ]{$ c } dx 


K] = o 

K] = o 

K] = o 
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[K 

] = J m[cosA s cosA a 

cos0 7 +sinA s 

cos 9] 


+ cos A s sinA s 

sin0; 

]{$?} dx 


K; 

1 We 

| = / m { 2 cos 9 1 

J Q 

cos A s 

cos A a + cos A s 

cos A a cti 


+ cosA* sinA a 

sin 9j 

+ sin 9 j sin A s 

]{*,} dx 

K] 

We 

= / — m[ sin 

•/ o 

cos A s 

sinA a + sinA s 

cos 6 1 + cosA s cosA 0 


+ sin0/ cosA a 

; sin A s 

]{$J dx 



ANTRAL LIBRARY 

1 *• T„ KAHPU* 

‘*" faA 12808 $ 



52 


6.2.7 Matrix [M 3 ]i 4x3 

[ M n\ = J Q m[c°sA 3 {-cosip k (-asindr + (x k + u k ) cosA 4 sinA a 
+t’fc sin A, sinA a + w k cosA a )} 

+ sinA s cosA a {—sin ip k cos dr (—a sin Or + (x k + u k ) cos A, sinA a 
+yjt sin A 4 sin A a + w k cos A a ) 

— sinV'fc sin Or (a cos 9 t — (x k + u k ) sin A, +v k cosA s )} 

+ sinA a sinA a (cos Tjj k (a cos Or — (x k + u k ) sin A, + v k cos A s ) 

n-1 

+ sin xp k cos 6 1 (ex + + 53 Ge)» + {x k + u k ) cos A, cos A a + ujt sin A, 

i=l 

cosAa -w k sinAa )}]{$ c } dx 



f 1 ' 

J m[cosA, {-sin^* (-asin0/ +(xjfe+Ufc) cos A, sinA a 
+v k sinA s sinA a +w k cosA a )} 

+ sinA 4 cosAa {costp k cos Or (— asin#/ + (x k + u k ) cosA 4 sinA a 

+v k sin A, sinA a + w k cosA fl ) - cos ip k 

sin 9 1 (a cos Or -(x k + u k ) sinA 4 +v k cos A, )} 

+ sin A s sinAa (sin^ (a cos Or — (x k + u k ) sin A, + v k cosA 4 ) 

n — 1 

- cos V'jfc cos Or (ex + e 2 + 0e)» + ( x k + u k ) cos A, cos A a + v k sin A s 

!=1 

cosA a - Wk sinAa )}]{$,} dx 



l 


h "Z* 

m[cosA, cos Or (ex + e 2 + 2 _,(/ e )i + ( x * + u*) cos A, cos A a +v k sin A, cosA a 

i=i 

-w k sin A^ ) — sin A 4 cos A a cos Or (a cti 


~(x k + u k ) sin A, +v k cosA 4 )+sinA 4 sinA a sin 0/ 

n-l 

(ex + e 2 + (l e )i + (x k + u k ) cos A, cos A a +v k sinA 4 cosAa 

»=i 
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-w k sin A a )]{<E» C } dx 

[<] = fo m [ cos A a {cos xl) k (a cos 9 1 — (x k + u k ) sin A, + v k cos A, ) 

+ sin cos 0/ (d + e 2 + ]T (l e ) { + (x fc + u k ) cos A, cos A a + v k sin A 4 cos A a 

i=X 

—w k sinA a ) + sinA 4 {cos0/ sinifik (— asin0/ 

+(** + «*) cosA, sinA a + v k sin A, sinA a +w k cos A. ) 

+ sin^>fc sin 0/ (acos0/ — (x^ + u*) sin A 4 + cosA 4 )}]{$ c } dx 

[M 2 3 2 ] = ^ m[cosA a {sin tp k (acos0/ - (x* + u fc ) sin A, 

71-1 

+v k cos A 4 ) - cos ipk cos 9 1 (ei + e 2 + (/,.), + (x fc + u k ) cos A, cos A a 

i = 1 

4-ujt sin A 4 cos A a - w k sin A a ) } — 

sin A, {cos ip k cos 0/ (-a sin 0/ +(x* + u Jfc ) cosA 4 sin A a 

+x* sin A, sin A a + ui fc cos A a ) 

+ cos V’jfc sin 0/ (a cos 0/ - (x k + u k ) sin A 4 + v k cos A 4 ) }] {$ c } dx 

[ 3 1 ^“*1 

■^23 J = m [cos A a {sin 0; (ei + e 2 4 - Y] (l e ), + (x k + u k ) cos A„ cos A„ 4 - v k sin A 4 cos 

0 i=l 

-uifc sinAa )} +sin A 4 {cos0/ (-a sin 9 r +(x k + u k ) cosA 4 sin A a 
+v k sinA 4 sinA a +vJ k cosA a )+cos0/ (a cos 9/ -(x k + u k ) sinA 4 
+v k cos A, )}]{$ c } dx 

K] = o 


K] = o 
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0 



+(x k + u k ) cos A, sin A a +v k sinA 4 sinA a +w k cosA 0 ) 

— simfi k sin dj (a cos 6j — (x k + u k ) sinA 4 + v k cos A s )} 

+ sin A 4 {cos tp k (—a sin#/ +{x k + u k ) 

cos A 4 sinA a + v k sinA 4 sinA a + w k cosA a )} 

+ cos A s sinA a {cos^* (a cos Oj - (x k + u k ) sin A, + v k cos A, ) 

n — 1 

+ sin tp k cos Qj ( e x + e 2 + (Qi + (x* + u k ) cos A s cos A a + v k 

*=1 

sin A 4 cosAa - w k sin A a )}]{$,} dx 


fie 

I 771 [cos A 5 cos A fl {cos^ cos Oj (— asin 9 j 

JO 

+(x k + u k ) cos A 4 sinAa +v k sinA 4 sinA a +w k cosA a ) 

+ sin 9 1 cos ijj k (a cos Gj - (x k + u k ) sin A s 4- v k cosA s )} 

+ sin A 4 {sin ip k (-a sin # ; + (x k + u k ) cos A, sinA a 
+v k sin A s sinA a + w k cosA a )} 

+ cos A 4 sinAa {sin^* (a cos G, - (x k + u k ) sinA 4 +v k cos A 4 ) 

n -1 

— cosifik cos 6] (ei + e 2 + ^(/ e )i + (x*: + u*) cos A s cosAa + v k sinA s 

1=1 

cosA a ~w k sinAa )}]{$,} dx 


flc 

/ to[cosA 4 cosA a {—cos (-a sin 6} 

J 0 


+(x k + u k ) cosA 4 sinA a +v k sinA 4 sinA Q + w k cosA a ) 

— sin#/ (a cos 0] - (x k + u k ) sinA 4 +v k cosA 4 )} 

n— 1 

— sin A, {cos 61 (ei + e 2 + (k)* + ( x k + u k ) cos A 4 cos A c 

i = 1 



55 


+ujt sinA g cosA a - sin A a )} + cos A 5 sin A, 

n-1 

{sin 9j (ej + + 'fT, (h)i + (%k + Uk) cos A s cos A a + Vk sin A, cos A a 

i= 1 

-Wk sin A q )}]{$ 9 } dx 
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6.2.8 Matrix [M 4 ] 14x3 

[m^i] = J [— mcosAj cosA a {a cos 9 r 

~((3 3 cos9[ - /^sin#/ + 1) (^(i e )i + x k ) - (ei +e 2 )}cos^ 

i=l 

n-1 

+m cos A, cos A a {(2(0* + ®* )sin^ fc }]{$ c } dx 
1=1 

[MfJ = f [—m cos A 4 cos A a a cos 
» 0 

n-1 

— (A cos ^ — PjsinOr + 1) (£(l e )i + x k ) — (ei + e 2 )} sini/>* 

i=l 

n-1 

-m cos A, cosAa {(£(i e )j + x k )cos i> k }]{$ c } dx 

1=1 

[Miz] = / [macosA, cos0, (2cos0j -sin0, ){$ c } {$ c } ]dx 

J 0 


[^21 ] = / [~^cosA a cos A, cos 9 1 ( 2 a sin ip k cos 9 r + (^(J e )i + a;* )cosij} k }{<£> c } ]dx 


K] = ' [msinA a cosA a cosA s (2acosV'it cos#/ — + £* )sin^>fc }{$ e } }dx 

Jo i=i 

fie n ~ 1 

[M 23 ] = [{— mcos A a cos 2 A* (^(Z e )i + x k ){Q P + (3 d cos 9[ - (3, sin 0/ - 0/ )}{$ c } da; 


K] = 0 


[-^32] — 0 



0 
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[ M «] - {-m els cos A a cos dj {a cos 8j - (ft cos 6j - ft sin Qj + 1 ) C£(h)i + x k ) 
~( e i + e 2 )}sinft. {$ 9 } ]rfx 

KJ = [m cZs cos 6j {acosfl, - (ftcosfl, - A sin Qj 4 - 1) (£&)< + ** ) 

1=1 

( e l + e 2)} COS Ipk {<ft} \dx 
r n~ 1 

[^43 J = m E{( e ! + e 2 ) + 2(53(Z e )i + zjt )} cos 6] cos A 5 cos 2 A a {$J ]dx 



59 


6.2.10 Nonlinear Vector [V JV£ ]i 4x1 

[^n L ] = J 0 —m[29j cosA s cosA a cos 0/ {$ c } {$ c } r {F} 

-K6>‘/ 2 cos 9 j cos A , + 0; cos A s 9,-9) ){$ c } {$ c } r {If} 

+9) cos 2 A a {«l> c } {$,} T {[/} + 20; cosA 0 {$ c } {$ 9 } r {U}]dx 

[ V 2 i L ] = l ~m[(9j -9j )cosA s {$ c } {$ c } r {F} 

J v 

-9j cos A, cos A a {$ c } {$ C } T {F} + 0' 7 2 cos 2 A a cosA s {$ c } {$ c } r {W}]dx 

[V£ L ] = o 

[C L ] = [ m[(—9] ) cos A s cos A a {$„} {$ c } r {F} 

J o 

-0i +0/ 2 )cos0 7 cosA s cosA a {$J {$ C } T {W} 

9 1 cos A s {$J {$ 9 } r {U}]dx 
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6.3 ELEMENT MATRICES ASSOCIATED WITH 
STRAIN ENERGY VARIATION 

The element matrices associated with the strain energy variation are derived by substi- 
tuting the assumed expressions for the displacement functions (Eq. 6.2), in the strain 
energy variation, 8Ui (Eq. 5.9), and carrying out the integration over the length of the 
beam element. The resulting variation of the strain energy can be written in the form; 

SU, = [*;£]{,} + {F e }} (6.8) 

Where 

{q} represents the vector of unknown degree of freedom. 

[K E ] is the linear stiffness matrix. { F E } is the nonlinear stiffness vector. 

The detailed expressions for the \K E ] and {F E } are as follows: 


6.3.1 Linear Stiffness Matrix [K E ] 

The linear stiffness matrix [ K E ] is given by the following sub-matrices: 

[K E ] = f lt [(W^os9g -E^^^G)mm T ]dx 

J 0 

[ a - e ] = j‘‘ [(ET^costfc -SJ^sin9 0 ){*"} {$"} T ]dx 
[A-,*] = f' [(Ejs^ro {t;j T + (rEAmm ] *■ 

KJ = jf Jc-raow} ] dx 
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[*»] = KY 

K &\ = jT [(E^cosSc ~m^sme c ){*"} {^} T ]dx 

*£] = to t + (rEimm ] * 

'k£\ = [(-E?Q{*;'} {«;'} r ] & 

*1] = KF 

Kg] = [/fg] T 

[*£] = £ P^M*;'} {i>;'} r + (tea0 2 ){v;} {*;} r 
+{*;} iKf + (t 2 ead ' 3 + gj){*;} {*;} T ]<& 

[*£] = £ [(-^ J D„){*i'> {*;'} T + (-r £ADi){t;j {4.;'} r ] & 

K,] = [* 5 F 

[*£] = l*£F 

Ks] = i^] T 
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*«] = r [(^){*i> TO T ] dr 


6.3.2 Simplification of Linear Stiffness Matrix [K E ] 

A simplified form of the linear stiffness matrix can be obtained by making the follow- 
ing assumptions: 


1 . The elastic axis is coincident with the modulus weighted centroid. 

2. The warping function is considered only in torsion. 

If the properties like El, EI m ,EI ^ etc., vary linearly along the element length, vari- 
ous quantities defined in the linear stiffness matrices can be expressed as 
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El « = 


Eh n + l z 


COS 

sin^c 


s/„ { = 


EI«, + p J ”() - Ehti \ x 


EJ m + f ^a - ^a ') i 


COS $g 

sin#G 


Ehn = 


£/ ((] + * 


+ 


El^m + l EI <"A x 


sin 0 G 
sin 0 g 


Elm = 


£/„ fl + ( ^ EI ** 


+ 


Elrm i + 


^ y 

Elrjrji ~ EJjjV2 


sin ^g 
cos 0 g 




E 4 A, + ( £ ‘ 4gl2 ~ £j4£ ’ 11 


) J 


cos 9 c 


EAD 2 1 H" 


/ EAD22 ~ EAD21 
V I . 


sin 0c 


£MZX 


£MD' U + 


EAD' 12 - EAD\ 


11 


x 


■ EAD , i+ ^A^EA^y 


cos 6g 
sin 6 g 


EAr) a 


fEAria2 - EArjgy 


EA Vn + { — ~‘ L ) * 

\ EAin+ (EA^m^ 


cos 9q 

sin0G 



64 


EXD 3 

EAD h 


ead 31 + (E4P**- E AD*1>) x ] 

EAD 5l + (EA Dt2-EAD^ x j 


EAD' 


EAD' 31 + 


^ EXE' 2 - E AD: 


51 


GJ 


EAD 0 


= GJ 01 + o' ) x ' 

^ Gc/22 — GJ21 ^ 


-2 


GJ U + ^J, 2 -GJ 11 ^ 


GJoi + 


EADn 1 H" 


( EAD 0 2 — EAD 1 


01 


x 


EX 


EX] + 


^EX 2 — EXi ^ 


Where the subscript 1 and 2 refer to the properties at the two end nodes. Substituting 
these expressions in the stiffness matrix leads to several integrals along the length of 
the element, which are given as 


12 6/ e -12 6/ f ' 

6/ e 4/2 —6/f 2/ t 2 

-12 -6/e 12 -6/e 

6/e 2/2 -6/e 4/2 . 


[ lt TO TO r ^ 

J 0 


1 

3/ e 


7 -8 1 
-8 16 8 
1 8 7 


f‘ TO TO 5 

Jo 


dx — 


ll 
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w) {*;r dx = 

C {*;') <f;r dx = 

/'{*;} m T dx = 

J‘‘ «} {fi'} T dx = 
/' {«;} {$;r dx = 
C *«} m T dx = 

Jo 

t m T d * = 

J 0 
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IS w <& 

IS *{*;'} {»;r dx 
IS *{«;'} {*;y dx 

IS *{*"} w r dx 

IS x{*;} {*;r dx 


1_ 

% 


1_ 

ll 


1 

3/ e 


1 

3/ e 


1 

3/ e 


4 

OO 

1 

4 ' 

0 

0 

0 

-4 

OO 

-4 

4/ e 

1 

<u 

<v 

00 

1 


8 

-16 

8 

-16 

32 

-16 

8 

-16 

8 


-2 

-8 

10 ‘ 


4 

16 

20 


-2 

-8 

10 _ 


3 

-12 

9 

1 

2/ e 

-4/ e 

2/e 

-3 

12 

-9 

7/ e 

-8/ e 

7/e . 

-2 

4 ■ 

-2 ‘ 


-8 

16 ■ 

-8 


10 

20 

10 



(6.9) 
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6.3.3 Nonlinear Stiffness Vector [.F-Ej 

The nonlinear stiffness vector [F E ] is given by the following sub-vectors: 

{Ff} = {[ (K)(^a) — (S' x sin 6^ )(v xx cos d G + w xx sin 0 G )] 

+[-(EI^, cos 9 G EI^ sin 6 g )(^)(«**) + (ET^cos9 g -E7^sm6 G ){<l>){w xx ) 
— - ( EArj a (u| +w 2 x ) + (S' x cos 8 G ) (~v x sin 0 G + w x cos 0 G ) 

+(M; C os0 g - M'(- sin 6 g )(<£) - l(£lC^)]{$''} }dx 

{Ff} = f { [ (K)(w x ) - (S' x cos e G ) (v xx cos d G + w xx sin B G )] 

J 0 

+[- (EI^ sin 6 G EI m cos 6 G )(<j>){v xx ) 4- {Eltsin 0 G + EI^ cos d G )(<j>)(w xx ) 
~ ^ (FAQ (v x + w x ) + (S' x sin 6 a ) (~v x sin d G + w x cos 6 G ) 

+(M' sin 6 G + M' cos0 G ){<j>) - ^{EAQ) (<%)}{$"} }dx 

{F?} = J {{(M'rt cos 0 G - M[ sin 9 G )(v xx ) + (M' sin 0 G + M[ cos 6 a )(to IX ){$,} 
[(r 0 )[(EAD^v xx - EAD[w xx m 

+\ead' 0 {vI 4- wl) + \ead'mI)} 

~(GJ 0 - GJ \)(4> 0 ) - (MMW 

-[^ EAD 0 {v\ + wl) + ]-EAD 4 {4> 2 x ) + (EAD 2 v xx - EAD : w xx (p]{^} }dx 

Z Z 

{if } = J {[(EAQv xx - EZifcM + \EA{vl + to*) + 1 -EAC 0 {4>1)){^} 

The underlined terms in {if}, {if} and {if} are associated with the axial strain 
at the elastic axis (u x + \v\ + \w\) of the blade. These nonlinear terms are modi- 
fied to a set of linear terms by substituting the axial stress by the axial inert ial force. 
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The description of the substitution approach [19, 29] is given in next section for clarity. 

6.3.4 Treatment of Nonlinear Terms Associated with Axial 
Strain at Elastic Axis 

The equation of motion corresponding to the axial degree of freedom can be written 
symbolically as [19,29]: 

m Q2 ^ *1* — Z u = 0 (6.10) 

The axial stress resultant (Eq. 5.10) can be rewritten as: 

V x = EA[u x + ^vl+ l -wl] + F ( 6 . 11 ) 

Where F contains all the additional higher order terms, which can be neglected. The 
simplified expressions for V x can be written as: 

% k EA\u, + 1 vl + (6.12) 

Neglecting all the higher order terms and also the time derivative terms, the dis- 
tributed inertia force Z u can be written as: 


where 


Z u = Cl[X^e)t +2jt] + Ci 

t=l 

c x = {l-(& cos 6j +(3 d sin 6 1 )}cos0/CosA s cosA a 

C 2 = {(ei + e 2 ) - a(p d sin 0/ + & cos 6j )} cos 6] cos A s cos A a 


(6.13) 


Integrating Eq. 6.10, 


Eo 

mCt 2 


i-L 


Z u dx 


■x k 


(6.14) 
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V x can be written as 


— /(£) — a]£ 2 + a2^ + 03 


(6.15) 


where 




Xk 


ai, 02 , a 3 are given as: 


_ rnQ 2 1 2 

01 - '£T [ ~2 ClZe] 


= 2p^[-c, £('=)< - cyi. 

i=l 


E n 


mO 2 1 1-L 

Z t=l i=l 


A T n-1 

%)i - 
t=i 1=1 


Combining Eqs. 6.15 with 6.12, the axial strain at the elastic axis can be written in 
terms of axial inertia force as 

k + \vl + \w 2 x ] = (6.16) 

Using Eq. 6.16, the underlined terms in {Ff}, {F 2 £ } and {Ff} are modified to a set 
of linear stiffness matrices, which are given as: 

[ 1 ‘ %v x mdx - i e f mmm T duv} = i kshv) 

Jo JO 
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r lc - 

/ V x w x {cj)'}dx 
Jo 

J o f x 4> x {4>' q }dx 


= [ /(OTOTO^H'} = [ Kg\{W] 

= 'w] '*/’ = [*£!W 


Combining these submatrices, a linear stiffness matrix [K E '] is obtained, which can 
be written as 


where 



[Kg’] [0] t0] (0] 

[0] [Kg] 10] ]0] 

[0] 10] [Kg] [0] 

[ 0 ] [ 0 ] [ 0 ] [ 0 ] 


(6.17) 


[K\i ) = [K§i ] = (°i + °2 + ° 3 )[-^-i] — (2a i + 02)^2] + (2ai ) [A3] 


[K§1] = (fli + a 2 + az)[Bi] - (2qi + a 2 )[B 2 ] + (2ai)[B 3 ] 


where 




’ 36 3 l e -36 3/ e 

3 l e 4/2 -3/e ~/2 

-36 -3/e 36 -3/e 

3/e -ll -3/e 4/ e 2 . 


’ 36 0 -36 6/ e 

0 6/2 0 - l 2 e 

-36 0 36 -6/ e 

6/e -ll -6/e 2/2 . 


(6.18) 


(6.19) 
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[ 4 *] = 


1 

' 24 

— 2l e 

-24 

5Z e 

-2 l e 

6/ e 2 

2l e 

-/ 2 

140 l e 

-24 

2l e 

24 

— 5/ e 


5 l e 

-ll 

-5 l e 

1/2 

3 e 


( 6 . 20 ) 




EACq 

3EAl e 


7-8 1 ‘ 

-8 16 -8 
1-8 7 


( 6 . 21 ) 



EAC 0 

6EAl e 


11 -12 1 
-12 16 -4 
1 -4 3 


( 6 . 22 ) 



EAC 0 
3 EAl e 


23 -26 3 

-26 32 -6 

3-6 3 


(6.23) 



Chapter 7 

RESULTS AND DISCUSSION 


The first step in any aeroelastic response and stability analysis is the evaluation of 
natural frequencies of the rotor blade. Using the inertial and structural model devel- 
oped in this study, a structural dynamic analysis was performed. It may be noted 
that the inertial and the structural operators given in Eqs. 6.6 and 6.8 respectively, 
are nonlinear. Since the structural dynamic analysis requires only linear terms, all 
the nonlinear terms are neglected. The corresponding linear equation for one beam 
finite element can be written as: 


+ = o (7.i) 

where [M]* represents the mass matrix of i th element, (given in Eq. 6.6. and in Sec. 
6.2.1). The stiffness matrix [K]i consists of three components. They are [K ci ] (given 
in Sec. 6.2.2), [K E ] (given in Sec. 6.3.1) and [K E '\ (given in Eq. 6.17). The element 
stiffness matrix is given as: 


[K\ i = [K CF ] + [K F ] + [K E '] (7.2) 

The element matrices are assembled to form the global finite element model for the 
rotor blade. For the tip element, the corresponding local to global transformation 
(Ref. 30) is performed to take into account the sweep and anhedral angles. The 
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transformation matrix [A c ] is given in Appendix A. Imposing the root boundary con- 
ditions, the corresponding rows and columns from the global model are eliminated. 
The resulting matrix equation can written as: 

[M]{q} + [K]{q} = 0 (7.3) 

Performing an eigen analysis, the natural frequencies of an undamped rotating blade 
in vacuum can be evaluated. 


7.1 VALIDATION 

In order to validate the finite element blade model developed in this study, the results 
of the present analysis are compared (for certain specific cases) with those available 
in the literature. The data shown in Table 1, correspond to a uniform, untwisted 
straight blade for which uncoupled natural frequencies are available in the literature 
(Refs. 24,26 and 27). The terminology soft-in-plane blade configuration indicates 
that the first non-dimensional rotating lag frequency is less than 1. 

Using the data given Table 1, the uncoupled natural frequencies of the rotating blade 
are calculated. A convergence study on the number of elements is also performed. 
The results of the convergence study along with the results of Refs. 24,26 and 27 are 
shown in Table 2. The results of the present analysis are in excellent agreement with 
those available in literature. The natural frequencies converge to a good accuracy 
with 10 elements. 

In all the subsequent calculations only 10 elements are considered. The effects vary- 
ing presweep (/?$), precone (/3 p ) and predroop (/?<*) on the natural frequencies of the 
straight rotor blade are shown in Table 3. The natural frequencies in flap, lag and 
torsional modes reduce with pre sweep angle, whereas they are not influenced by pre- 
cone and predroop angles. 
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7 .2 EFFECTS OF TIP SWEEP AND ANHEDRAL 

The influence of tip sweep and anhedral angles on the natural frequencies and mode 
shapes of the rotor blade are analysed systematically. The blade data is given in 
Table 1. The blade is divided into 10 finite elements. The swept tip portion of 
the blade is taken as one element having a length 0.11. The sweep and anhedral 
angles are varied from —30° to +30° in steps of 5 degrees. For each case the coupled 
natural frequencies and mode shapes are evaluated. In total, there are 169 cases. The 
frequencies corresponding to all these cases are given in Appendix B, for reference 
and for future validation purposes. 

The variation of natural frequencies with sweep angles are shown in 3-D plots in 
Figs. 12 to 17. Figure 12 shows the variation of first lag frequency with sweep 
angles. The results show that depending on the value of anhedral angle, the first lag 
frequency decreases monotonically with sweep angle ( for A a = +30° ) or decreases 
and increases with sweep angles (for A a = -30°). The first lag frequency varies from 
0.8898 (A a = 30°, A s = —30°) to 0.6892 (A a = 25°, A s = 30°). A similar nature is also 
observed for second lag frequency, as shown in Fig. 13. 

The variation of first flap frequency with sweep angles is shown in Fig. 14. It is 
interesting to note that depending on the value of anhedral angle, the flap frequency 
increases monotonically (for A 0 = +30°) or decreases monotonically (for A fl = —30°) 
with sweep angles. The influence of sweep angle is more pronounced for large values 
of anhedral angles. The first flap frequency varies from 1.2036 (A a = A s = 30°) to 
1.1003 (A a = 30°, A 5 = -30°). Similar behaviour is also seen for second flap fre- 
quency, shown in Fig. 15. The third flap frequency exhibits a different phenomenon 
ith variation in sweep and anhedral angles (Fig. 16). The reason may be due to 
complex modal coupling with lag and torsional modes. The torsional frequency 
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seems to be less influenced by sweep and anhedral angles as shown in Fig. 17. 

It is found from the results that anhedral angle alone introduces lag-torsion coupling 
and sweep angle introduces flap-torsion coupling. The influence tip anhedral angle on 
the mode shape in first lag is shown in Fig. 18. The torsional mode contribution is 
increasing with increase in the value of A a . The flap motion has negligible contribu- 
tion in first lag mode. 

The influence of tip sweep on first flap mode is shown in Fig. 19. Even in this case, 
the participation of torsional motion increases with increase in A s . The participation 
of lag motion is negligible. 

The combined effect of both anhedral and sweep angles on the mode shapes are shown 
in Figs. 20-25. The combination of sweep and anhedral angles introduce significant 
lag-flap-torsion coupling in first lag mode (Fig. 20), whereas the flap contribution is 
negligible in second lag mode (Fig. 21). The contribution of lag mode is negligible in 
all the flap modes (Figs. 22-24) even with influence of A a . Torsional mode shape is 
unaffected by sweep and anhedral angles as shown in Fig. 25. 

To quantify the lag-torsion and flap-torsion coupling in various modes, the value of 
torsional deformation at 0.75 1 ($ 0 . 75 ) is taken as the reference parameter. (Since 
this cross-section is treated as a typical section for aerodynamic load calculations in a 
rotor blade). The variation $ 0.75 as a function of sweep angles is shown in Figs. 26-30 
for the lag and flap modes. 

In lag modes, lag-torsion coupling is significantly influenced by anhedral angle com- 
pared to sweep angle as shown in Figs. 26 and 27. In first flap mode, the sweep angle 
has a significant influence on flap-torsion coupling compared to anhedral angle. In 
second flap mode, both anhedral and sweep angles influence the flap-torsion coupling. 
A similar variation is observed in third flap mode but of a much smaller magnitude. 



Chapter 8 

CONCLUDING REMARKS 


In this thesis, finite element model of a rotor blade with swept tip has been developed 
for the dynamic analysis of the blade. The model includes all the geometric complex- 
ities of the hub and hub motion. 

The equations of motion have been derived using Hamilton’s principle . The general 
equations are then specialised for the structural dynamic analysis of the rotor blade. 
The formulation was validated by comparing the results of the present analysis (for a 
straight blade) with those available in literature. 

The influence of tip sweep on natural frequencies and mode shapes has been analysed 
systematically. Anhedral angle introduces significant lag-torsion coupling and sweep 
angle introduces flap-torsion coupling. 
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Table 1: Input data for soft-in-plane blade 

Soft-in-plane blade data 

Trn^JmP = 0.004 
Im^/ml 2 = 0.000 
6g — 0 
m = 1 
& = 0 

Pd = 0 

P P = 0 

Or = 0 

GJ/mtt 2 l 4 = 0.001473 
EI/mQ 2 l 2 = 20.0 
e 2 = 0.0 
e 2 = 0.0 
a = 0.0 

EACJEA = 0.00021036 
EI a /mQ, 2 l 4 = 0.0301 
EJrJmSW = 0.0106 


Offsets of mass centre and aerodynamic centre from elastic centre are zero. 
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Table 2: Natural frequencies of soft-in-plane blade 


Mode 

Uncoupled Natural Frequencies 

No. of elements 

“Ref. [26] 

Ref. [27] 

Ref. [24] 

N=2 

N=5 

N=10 

N=20 

1st Lag 

0.73519 

0.73128 

0.73114 

0.73113 

0.7326 


0.732 

2nd lag 

4.49322 

4.45488 

4.45317 

4.45306 

4.4563 

4.4825 

- 

1st Flap 

1.13260 

1.12697 

1.12578 

1.12513 

1.1247 

1.1245 

1.125 

2nd Flap 

3.45728 

3.44530 

3.44239 

3.42668 

3.4089 


_ 

3rd Flap 

7.79626 

7.77756 

7.77462 

7.71541 

7.6376 

7.6218 

- 

1st Torsion 

3.26406 

3.26330 

3.26328 

3.26328 



3.263 

1st Axial 

6.94918 

6.94035 

6.93925 

M5. 93898 

6.9533 

- 

- 


Table 3: Influence of precone, predroop and preswept angles on natural frequencies 


Mode 


& = 2 

(3 S =3 

A> = 2 

= 2 




(3 P = 0 

Ps = 0 





13d =0 

(3 d = 0 


1st Lag 


0.7034 

0.6891 

0.7311 

imu 

2nd Lag 

tfwal 

4.4280 

4.4154 

4.4531 


1st Flap 

1.1257 

1.1085 

1.0998 


1.1257 

2st Flap 

3.4423 

3.4088 

3.3920 

3.4423 

3.4423 

3st Flap 

7.7746 

7.7334 

7.7127 

7.7746 

7.7746 

1st Torsion 

3.2632 

3.2602 

3.2587 

3.2632 1 

3.2632 

1st Axial 

6.9392 

6.9392 

6.9392 

6.9392 

6.9392 
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Figure 12: Non-dimensional first lag natural frequency as a function of tip sweep angle 


second lag frequency 
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Figure 13: Non-dimensional second lag natural frequency as a function of tip sweep 
angle 


first flap frequency 



Figure 14: Non-dimensional first flap natural frequency as a function of tip sweep 
angle 


second flap frequency 
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Figure 15: Non-dimensional second flap natural frequency as a function of tip sweep 
angle 


third flap frequency 


92 



Figure 16 : Non-dimensional third flap natural frequency as a function of tip sweep 
angle 


torsion frequency 
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Figure 17: Non-dimensional torsion natural frequency as a function of tip sweep angle 
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Figure 26: Lag-torsion coupling Vs tip sweep angle in first lag mode 


0.75 
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Figure 27. Lag-torsion coupling Vs tip sweep angle in second lag mode 



0.75 
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Figure 28: Flap-torsion coupling Vs tip sweep angle in first flap mode 


0.75 
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Figure 29: Flap-torsion coupling Vs tip sweep angle in second flap mode 




References 


[1] Houbolt, J.C., and Brooks, G.W., ’’Differential Equations of Motion for Com- 
bined Flapwise Bending, Chordwise Bending and Torsion of Twisted Nonuni- 
form Rotor Blades,” NACA Report 1346, 1958. 

[2] Hodges, D.H. and Dowell, E.H., ’’Nonlinear Equations of Motion for the Elastic 
Bending and Torsion of Twisted Nonuniform Rotor Blades,” NASA TN D-7818, 
December 1974. 

[3] Rosen, A., and Friedmann, P.P., ’’Nonlinear Equations of Equilibrium for Elastic 
Helicopter or Wind Turbine Blades Undergoing Moderate Deflection,” NASA 
CR-1 59478, October 1978. 

[4] Kaza, K.R. and Kvaternik, R.G., ’’Nonlinear Aeroelastic Equations for Com- 
bined Flapwise Bending, Chordwise Bending, Torsion and Extension of Twisted 
Nonuniform Rotor Blades in Forward Flight,” NASA TM-74059, 1977. 

[5] Crespo DaSilva, M.R.M. and Hodges, D.H., "Nonlinear Flexure and Torsion of 
Rotating Beams with Application to Helicopter Blades-I. Formulation,” Vertica, 
Vol. 10, 1986, pp. 151-169. 

[6] Crespo DaSilva, M.R.M. and Hodges, D.H.. ’’Nonlinear Flexure and Torsion of 
Rotating Beams with Application to Helicopter Blades-II. Results for Hover," 
Vertica, Vol. 10, 1986, pp. 171-186. 


107 



108 


[7] Rosen, A. and Friedmann, P.R, ’’The Nonlinear Behavior of Elastic Slender 
Straight Beams Undergoing Small Strains and Moderate Rotations,” Journal of 
Applied Mechanics, Vol. 46, March 1979, pp. 161-168. 

[8] Hodges, D.H. and Ormiston, R.A., ’’Stability of Elastic Bending and Torsion of 
Uniform Cantilever Rotor Blades in Hover With Variable Structural Coupling,” 
NASA TN D-8192, 1976. 

[9] Shamie, J. and Friedmann, P.P., ’’Effect of Moderate Deflections on the Aeroe- 
lastic Stability of a Rotor Blade in Forward Flight,” Proc. 3rd European Rotor- 
craft and Powered Lift Aircraft Forum, Aix-en-Provence, 1977, pp. 24. 1-24. 37. 

[10] Sivaneri, N.T. and Chopra, I., ’’Dynamic Stability of a Rotor Blade Using Finite 
Element Analysis,” AIAA Journal, Vol. 20, No. 5, 1982, pp. 716-723. 

[11] Celi, R. and Friedmann, P.P., ’’Rotor Blade Aeroelasticity in Forward Flight 
with an Implicit Aerodynamic Formulation,” AIAA Journal, Vol. 26, No. 2, 
1988, pp. 1425-1433. 

[12] Bauchau, O.A., ”A Beam Theory for Anistropic Materials,” Journal of Applied 
Mechanics, Vol. 52, June 1985, pp. 416-422. 

[13] Bauchau, O.A. and Hong, C.H., ’’Finite Element Approach to Rotor Blade 
Modelling,” Journal of American Helicopter Society, Vol. 32, No. 1, Jan. 1987, 
pp. 60-67. 

[14] Bauchau, O.A., Coffenberrv, B.S. and Rehfield. LAV, ’’Composite Box Beam 
Analysis: Theory and Experiments,” Journal of Reinforced Plastic and Com- 
posites, Vol. 6, 1987, pp. 25-35. 



109 


[15] Bauchau, O.A. and Hong, C.H., ’’Large Displacement Analysis of Naturally 
Curved and Twisted Composite Beams,” AIAA Journal, Vol. 25, No. 11, Nov. 
1987, pp. 1469-1475. 

[16] Bauchau, O.A. and Hong, C.H., ’’Nonlinear Composite Beam Theory,” Journal 
of Applied Mechanics, Vol. 55, March 1988, pp. 156-163. 

[17] Kosmatka, J.B. and Friedmann, P.P., ’’Vibration Analysis of Composite Tur- 
bopropellers Using a Nonlinear Beam- Type Finite Element Approach,” AIAA 
Journal, Vol. 27, No. 11, 1989, pp. 1606-1614. 

[18] Hodges, D.H., ’’Review of Composite Rotor Blade Modelling,” AIAA Journal, 
Vol. 28, No. 3, 1990, pp. 561-564. 

[19] Yaun, K.A., Venkatesan, C., and Friedmann, P.P., ’’Structural Dynamic Model 
of Composite Rotor Blades Undergoing Moderate Deflections,” in Recent Ad- 
vances in Structural Dynamic Modelling of Composite Rotor Blades and Thick 
Composites, ASME AD- Vol. 30, pp. 127-155, ASME Winter Annual Meeting, 
November 1992. 

[20] Smith, E.C. and Chopra. I., ’’Formulation and Evaluation of an Analytical 
Model for Composite Box Beams,” Journal of the American Helicopter Soci- 
ety, Vol. 36, No. 3, 1991, pp. 23-35. 

[21] Celi, R. and Friedmann, P.P., ” Aeroelastic Modeling of Swept Tip Rotor Blades 
Using Finite Elements,” Journal of the American Helicopter Society, Vol. 33, 
No. 2, April 1988, pp. 43-52. 

[22] Celi, R. and Friedmann, P.R, ’’Structural Optimization with Aeroelastic Con- 
straints of Rotor Blades with Straight and Swept tips,” AIAA Journal, Vol. 28, 
No. 5, 1990, pp. 928-936. 



110 


[23] Kim, K.C., and Chopra, I., ’’Aeroelastic Analysis of Swept, Anhedral, and Ta- 
pered Tip Rotor Blades,” Journal of the American Helicopte Society, Vol. 37, 
No. 1, 1992, pp. 15-30. 

[24] Yuan, K.A., Friedmann, P.P. and Venkatesan, C., ”A New Aeroelastic Model 
for Composite Rotor Blades with Straight and Swept Tips,” AIAA Paper No. 
92-2259, Proc. 33rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dy- 
namics and Materials Cong., Dallas, TX, April 1992, pp. 1371-1390. 

[25] Yen, J.G., ’’Effects of Blade Tip on Dynamics, Cost, Weight, Aerodynamic 
Performance and Aeroelastic Response,” Journal of the American Helicopter 
Society, Vol. 39, No. 4, Oct. 1994, pp. 37-45. 

[26] Gupta, P.K., ’’Structural Dynamics of Helicopter Rotor having Precone - 
Presweep - Predroop- Pretwist and Torque offset including Hub Motions”, 
M.Tech Thesis, Dept, of Aerospace Engineering. I.I.T Kanpur, July 1996. 

[27] Pohit, G., ’’Dynamics of a Bearingless Helicopter Rotor Blade with a Non-linear 
Elastomatic Constraint,” Ph.D Thesis, Dept, of Mechanical Engineering. I.I.T 
Kanpur, Feb. 1999. 

[28] Yuan, K.A and Friedmann, P:P., ’’Aeroelastic and Structural Optimization of 
Composite Helicopter Rotor Blades With Swept Tips,” NASA Contractor Re- 
port 4665, May 1995. 

[29] Venkatesan. C., ’’Treatment of Axial Mode in Structural Dynamics and Aeroe- 
lastic Analysis of Rotor Blade -A Review,” Computational Structural Mechan- 
ics, Proceedings of the NASAS-94, I.I.T. Kharagpur, Dec. 8-9, 1994, pp. 361-372. 

[30] Panda. B., ’’Technical Note: Assembly of Moderate-Rotation Finite Elements 
Used in Helicopter Rotor Dynamics,” Journal of the American Helicopter Soci- 
ety, Vol. 32, No. 4, 1987, pp. 63-69. 



Appendix A 

Local to global coordinate 
transformation for tip element 


The linear transformation matrix to modify the local coordinate system of the tip 
element to global system is given by the expression. 


{*<?} = [A L ]tf} 

where {<?*} refer to the nodal degrees of freedom of the tip element in element coor- 
dinate system and {<7 G } corresponds to the global coordinate system. The non zero 
elements of [A L ] are given below. 



A L (1,1) 

= 

cos A 4 

A L (1, 12) 

= 

sin A 4 cosA a 

A L (2, 9) 

= 

— sin A 0 

A l (3,7) 

= 

sinA s sinA a 

A x (4,4) 

= 

cos A a 

A 1 ^) 

= 

cos A a 

A l (6,2) 

= 

-sinAa sinA £ 

A L (6, 9) 

= 

-sinA s cos A, 

A l (7, 14) 

= 

- sin A a 

A l (8, 6) 

= 

cos A s 

A l (9, 2) 

= 

cos A s sinA a 

A l (9, 9) 

= 

cosA s cosA a 

A^CIO, 4) 

= 

4 cos A s sinA a 

A L (10, 8) 

=r 

\ sin A s 

A L (11,4) 

= 

cos A s sinA 0 

A L (11, 11) 

= 

cos A s cosA a 

A L (l2,h) 

= 

cos A s sinAa 

A i (13,l) 

= 

-\ sin A s 

A z, (13, 5) 

= 

4 cos A s sin A a 

A i (13, 13) 

= 

cosA s cosAq 

A L (14,7) 


cos A 5 sinAa 


A L (1, 5) 

= 

sin A 4 sinAa 

A l (2, 2) 

= 

cosA a 

A i (3, 3) 

= 

cosAa 

A l (3, 14) 

= 

sinA 4 cosAa 

A z '(4, 11) 

= 

- sin A a 

A i (5, 12) 

= 

— sin A a 

A i (6,6) 

= 

cosAa 

A l (7, 7) 

= 

cos A a 

A l (8,4) 

= 

— sin A s sinA< 

A l (8,11) 

= 

— sinAa cos A, 

A l (9,6) 

= 

sin A s 

A £ (10,2) 

= 

| cos A s sinA 0 

A L (10, 6) 

= 

| sin A s 

A i (10, 10) 

= 

cosAa cosAa 

A L (11, 8) 

= 

sinAa 

A l (12, 1) 

= 

— sin A s 

A l (12, 12) 

= 

COS A s cosAa 

A L (13,3) 

= 

-4 sin A s 

A^nj) 

= 

| cos A s sin A fl 

A L (14,3) 

= 

- sin A a 

A l (14, 14) 

Z= 

cos A* cosAa 



Appendix B 

Natural frequencies and coupling 
measures of rotor blade with swept 
tips 
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Table B.l: Natural frequencies 


A a 

deg 

Aj 

deg 

1 lag 

1 flap 

torsion 

2 flap 

2 lag 

axial 

3 flap 

-30. 

-30. 

.8010 

1.2015 

3.2537 

3.7440 

4.7577 

7.1196 

8.2532 

-30. 

-25. 

.7820 

1.1994 

3.2542 

3.7163 

4.6724 

7.0912 

8.1994 

-30. 

-20. 

.7658 

1.1961 

3.2547 

3.6880 

4.6025 

7.0745 

8.1470 

-30. 

-15. 

.7522 

1.1915 

3.2551 

3.6594 

4.5463 

7.0695 

8.0965 

-30. 

-10. 

.7415 

1.1858 

3.2554 

3.6307 

4.5024 

7.0762 

8.0480 

-30. 

-5. 

.7336 

1.1789 

3.2555 

3.6022 

4.4694 

7.0937 

8.0019 

-30. 

0. 

.7284 

1.1709 

3.2556 

3.5740 

4.4463 

7.1211 

7.9581 

-30. 

5. 

.7259 

1.1619 

3.2555 

3.5463 

4.4325 

7.1572 

7.9168 

-30. 

10. 

.7259 

1.1521 

3.2554 

3.5193 

4.4272 

7.2005 

7.8779 

-30. 

15. 

.7283 

1.1415 

3.2551 

3.4929 

4.4300 

7.2492 

7.8413 

-30. 

20. 

.7330 

1.1304 

3.2547 

3.4675 

4.4406 

7.3015 

7.8069 

-30. 

25. 

.7397 

1.1188 

3.2542 

3.4431 

4.4586 

7.3556 

7.7745 

-30. 

30. 

.7481 

1.1071 

3.2537 

3.4197 

4.4838 

7.4099 

7.7440 

-25. 

-30. 

.8059 

1.1841 

3.2564 

3.6600 

4.8149 

7.1008 

8.1052 

-25. 

-25. 

.7866 

1.1819 

3.2568 

3.6399 

4.7156 

7.0673 

8.0693 

-25. 

-20. 

.7699 

1.1788 

3.2572 

3.6192 

4.6345 

7.0449 

8.0336 

-25. 

-15. 

.7558 

1.1748 

3.2574 

3.5979 

4.5691 

7.0341 

7.9983 

-25. 

-10. 

.7444 

1.1699 

3.2576 

3.5763 

4.5173 

7.0347 

7.9638 

-25. 

-5. 

.7355 

1.1641 

3.2578 

3.5545 

4.4775 

7.0463 

7.9302 

-25. 

0. 

.7292 

1.1575 

3.2578 

3.5326 

4.4484 

7.0680 

7.8977 

-25. 

5. 

.7254 

1.1502 

3.2578 

3.5109 

4.4288 

7.0990 

7.8665 

-25. 

10. 

.7238 

1.1422 

3.2576 

3.4894 

4.4179 

7.1379 

7.8365 

-25. 

15. 

.7244 

1.1337 

3.2574 

3.4683 

4.4150 

7.1833 

7.8078 

-25. 

20. 

.7270 

1.1248 

3.2572 

3.4476 

4.4196 

7.2335 

7.7804 

-25. 

25. 

.7315 

1.1156 

3.2568 

3.4275 

4.4312 

7.2870 

7.7543 

-25. 

30. 

.7375 

1.1063 

3.2564 

3.4080 

4.4495 

7.3421 

7.7295 
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- 20 . 

- 30 . 

.8113 

1.1686 

3.2588 

3.5913 

4.8743 

7.0884 

7.9936 

- 20 . 

- 25 . 

.7915 

1.1665 

3.2590 

3.5775 

4.7598 

7.0502 

7.9704 

- 20 . 

- 20 . 

.7741 

1.1637 

3.2593 

3.5629 

4.6666 

7.0226 

7.9466 

- 20 . 

- 15 . 

.7594 

1.1603 

3.2595 

3.5477 

4.5915 

7.0062 

7.9227 

- 20 . 

- 10 . 

.7471 

1.1562 

3.2596 

3.5320 

4.5318 

7.0010 

7.8987 

- 20 . 

- 5 . 

.7373 

1.1515 

3.2597 

3.5159 

4.4853 

7.0068 

7.8748 

- 20 . 

0. 

.7299 

1.1463 

3.2597 

3.4995 

4.4501 

7.0230 

7.8512 

- 20 . 

5 . 

.7248 

1.1406 

3.2597 

3.4829 

4.4250 

7.0488 

7.8281 

- 20 . 

10 . 

.7217 

1.1344 

3.2596 

3.4663 

4.4087 

7.0832 

7.8054 

- 20 . 

15 . 

.7207 

1.1278 

3.2595 

3.4498 

4.4006 

7.1251 

7.7834 

- 20 . 

20 . 

.7215 

1.1210 

3.2593 

3.4334 

4.3998 

7.1730 

7.7620 

- 20 . 

25 . 

.7240 

1.1140 

3.2590 

3.4172 

4.4058 

7.2254 

7.7413 

- 20 . 

30 . 

.7279 

1.1068 

3.2588 

3.4014 

4.4180 

7.2809 

7.7213 

- 15 . 

- 30 . 

.8172 

1.1549 

3.2607 

3.5360 

4.9354 

7.0832 

7.9092 

- 15 . 

- 25 . 

.7966 

1.1530 

3.2609 

3.5273 

4.8044 

7.0408 

7.8952 

- 15 . 

- 20 . 

.7785 

1.1507 

3.2610 

3.5179 

4.6986 

7.0084 

7.8804 

- 15 . 

- 15 . 

.7629 

1.1480 

3.2611 

3.5078 

4.6135 

6.9868 

7.8650 

- 15 . 

- 10 . 

.7497 

1.1449 

3.2612 

3.4971 

4.5457 

6.9761 

7.8492 

- 15 . 

- 5 . 

.7390 

1.1413 

3.2612 

3.4858 

4.4924 

6.9762 

7.8331 

- 15 . 

0. 

.7304 

1.1374 

3.2612 

3.4742 

4.4514 

6.9869 

7.8168 

- 15 . 

5 . 

.7241 

1.1332 

3.2612 

3.4622 

4.4210 

7.0076 

7.8003 

- 15 . 

10 . 

.7197 

1.1287 

3.2612 

3.4500 

4.3999 

7.0375 

7.7839 

- 15 . 

15 . 

.7172 

1.1239 

3.2611 

3.4376 

4.3870 

7.0756 

7.7676 

- 15 . 

20 . 

.7164 

1.1190 

3.2610 

3.4251 

4.3813 

7.1209 

7.7515 

- 15 . 

25 . 

.7171 

1.1139 

3.2609 

3.4125 

4.3823 

7.1719 

7.7356 

- 15 . 

30 . 

.7193 

1.1088 

3.2607 

3.4000 

4.3892 

7.2273 

7.7199 

- 10 . 

- 30 . 

.8237 

1.1430 

3.2621 

3.4921 

4.9976 

7.0858 

7.8455 

- 10 . 

- 25 . 

.8020 

1.1416 

3.2622 

3.4878 

4.8491 

7.0398 

7.8385 

- 10 . 

- 20 . 

.7829 

1.1400 

3.2623 

3.4827 

4.7300 

7.0031 

7.8307 

- 10 . 

- 15 . 

.7663 

1.1381 

3.2623 

3.4770 

4.6348 

6.9766 

7.8221 

- 10 . 

* 

J— ‘ 

O 

.7522 

1.1360 

3.2623 

3.4707 

4.5589 

6.9606 

7.8130 

- 10 . 

- 5 . 

.7404 

1.1336 

3.2624 

3.4638 

4.4990 

6.9553 

7.8033 

- 10 . 

0. 

.7308 

1.1310 

3.2624 

3.4564 

4.4524 

6.9606 

7.7931 

- 10 . 

5 . 

.7233 

1.1282 

3.2624 

3.4486 

4.4171 

6.9762 

7.7825 

- 10 . 

10 . 

.7178 

1.1252 

3.2623 

3.4403 

4.3914 

7.0014 

7.7715 

- 10 . 

15 . 

.7140 

1.1221 

3.2623 

3.4318 

4.3741 

7.0356 

7.7604 

- 10 . 

20 . 

.7118 

1.1189 

3.2623 

3.4229 

4.3642 

7.0779 

7.7490 

- 10 . 

25 . 

.7111 

1.1156 

3.2622 

3.4138 

4.3608 

7.1270 

7.7375 

- 10 . 

30 . 

.7117 

1.1123 

3.2621 

3.4045 

4.3632 

7.1817 

7.7260 
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- 5 . 

- 30 . 

.8308 

- 5 . 

- 25 . 

.8077 

- 5 . 

- 20 . 

.7874 

- 5 . 

- 15 . 

.7697 

- 5 . 

- 10 . 

.7545 

- 5 . 

- 5 . 

.7417 

- 5 . 

0 . 

.7311 

- 5 . 

5 . 

.7225 

- 5 . 

10 . 

.7159 

- 5 . 

15 . 

.7110 

- 5 . 

20 . 

.7078 

- 5 . 

25 . 

.7059 

- 5 . 

30 . 

.7053 

0 . 

- 30 . 

.8384 

0 . 

- 25 . 

.8137 

0 . 

- 20 . 

.7919 

0 . 

- 15 . 

.7729 

0 . 

- 10 . 

.7566 

0 . 

- 5 . 

.7427 

0 . 

0 . 

.7311 

0 . 

5 . 

.7217 

0 . 

10 . 

.7142 

0 . 

15 . 

.7084 

0 . 

20 . 

.7042 

0 . 

25 . 

.7015 

0 . 

30 . 

.7001 

5 . 

- 30 . 

.8466 

5 . 

- 25 . 

.8199 

5 . 

- 20 . 

.7965 

5 . 

- 15 . 

.7761 

5 . 

- 10 . 

.7585 

5 . 

- 5 . 

.7436 

5 . 

0 . 

.7311 

5 . 

5 . 

.7208 

5 . 

10 . 

.7125 

5 . 

15 . 

.7061 

5 . 

20 . 

.7013 

5 . 

25 . 

.6980 

5 . 

30 . 

.6959 


1.1329 3.2630 

1.1322 3.2630 
1.1315 3.2630 
1.1306 3.2630 
1.1295 3.2630 
1.1284 3.2631 
1.1271 3.2631 

1.1257 3.2631 

1.1242 3.2630 
1.1226 3.2630 
1.1210 3.2630 
1.1192 3.2630 
1.1174 3.2630 

1.1243 3.2633 
1.1248 3.2633 
1.1251 3.2633 
1.1254 3.2633 

1.1256 3.2633 

1.1257 3.2633 

1.1258 3.2633 

1.1257 3.2633 

1.1256 3.2633 
1.1254 3.2633 
1.1251 3.2633 
1.1248 3.2633 

1.1243 3.2633 

1.1172 3.2630 
1.1191 3.2630 
1.1209 3.2630 
1.1226 3.2630 
1.1242 3.2630 

1.1257 3.2631 

1.1271 3.2631 

1.1284 3.2631 
1.1295 3.2630 
1.1306 3.2630 
1.1315 3.2630 

1.1323 3.2630 
1.1329 3.2630 


3.4581 5.0599 
3.4576 4.8932 
3.4564 4.7606 
3.4547 4.6550 
3.4523 4.5712 
3.4494 4.5049 
3.4459 4.4530 
3.4419 4.4130 
3.4374 4.3832 
3.4325 4.3620 
3.4271 4.3483 
3.4213 4.3412 
3.4152 4.3397 

3.4328 5.1216 
3.4357 4.9361 
3.4381 4.7898 
3.4400 4.6741 
3.4413 4.5825 
3.4421 4.5101 
3.4424 4.4532 
3.4421 4.4090 
3.4413 4.3754 
3.4400 4.3508 
3.4381 4.3338 
3.4357 4.3234 
3.4328 4.3186 

3.4152 5.1815 
3.4213 4.9771 
3.4271 4.8173 
3.4325 4.6917 
3.4374 4.5927 
3.4419 4.5145 
3.4459 4.4530 
3.4494 4.4049 
3.4523 4.3680 
3.4547 4.3403 
3.4564 4.3205 
3.4576 4.3073 
3.4581 4.2999 


7.0968 

7.7980 

7.0476 

7.7966 

7.0071 

7.7945 

6.9761 

7.7916 

6.9551 

7.7881 

6.9446 

7.7840 

6.9446 

7.7792 

6.9550 

7.7739 

6.9755 

7.7680 

7.0056 

7.7616 

7.0445 

7.7549 

7.0912 

7.7477 

7.1447 

7.7402 

7.1163 

7.7636 

7.0647 

7.7669 

7.0209 

7.7697 

6.9857 

7.7718 

6.9601 

7.7734 

6.9445 

7.7743 

6.9393 

7.7746 

6.9445 

7.7743 

6.9601 

7.7734 

6.9857 

7.7718 

7.0209 

7.7697 

7.0647 

7.7669 

7.1163 

7.7636 

7.1447 

7.7401 

7.0912 

7.7477 

7.0445 

7.7549 

7.0056 

7.7616 

6.9755 

7.7680 

6.9550 

7.7739 

6.9446 

7.7792 

6.9446 

7.7840 

6.9551 

7.7881 

6.9761 

7.7916 

7.0071 

7.7945 

7.0476 

7.7966 

7.0968 

7.7980 
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10. 

-30. 

.8551 

1.1115 

3.2621 

3.4045 

5.2383 

7.1817 

7.7260 

10. 

-25. 

.8262 

1.1152 

3.2622 

3.4138 

5.0155 

7.1270 

7.7375 

10. 

-20. 

.8009 

1.1187 

3.2623 

3.4229 

4.8425 

7.0779 

7.7490 

10. 

-15. 

.7790 

1.1221 

3.2623 

3.4318 

4.7075 

7.0356 

7.7604 

10. 

-10. 

.7602 

1.1252 

3.2623 

3.4403 

4.6015 

7.0014 

7.7715 

10. 

-5. 

.7442 

1.1282 

3.2624 

3.4486 

4.5181 

6.9762 

7.7825 

10. 

0 . 

.7308 

1.1310 

3.2624 

3.4564 

4.4524 

6.9606 

7.7931 

10. 

5. 

.7198 

1.1336 

3.2624 

3.4638 

4.4009 

6.9553 

7.8033 

10. 

10. 

.7110 

1.1360 

3.2623 

3.4707 

4.3611 

6.9606 

7.8130 

10. 

15. 

.7040 

1.1382 

3.2623 

3.4770 

4.3308 

6.9766 

7.8221 

10. 

20. 

.6988 

1.1401 

3.2623 

3.4827 

4.3085 

7.0031 

7.8306 

10. 

25. 

.6952 

1.1419 

3.2622 

3.4878 

4.2931 

7.0398 

7.8384 

10. 

30. 

.6929 

1.1434 

3.2621 

3.4921 

4.2833 

7.0858 

7.8454 

15. 

-30. 

.8639 

1.1071 

3.2607 

3.4001 

5.2908 

7.2272 

7.7201 

15. 

-25. 

.8325 

1.1130 

3.2609 

3.4126 

5.0503 

7.1719 

7.7357 

15. 

-20. 

.8052 

1.1185 

3.2610 

3.4251 

4.8651 

7.1209 

7.7515 

15. 

-15. 

.7817 

1.1237 

3.2611 

3.4376 

4.7214 

7.0756 

7.7676 

15. 

-10. 

.7616 

1.1286 

3.2612 

3.4500 

4.6091 

7.0375 

7.7839 

15. 

-5. 

.7446 

1.1332 

3.2612 

3.4622 

4.5208 

7.0076 

7.8003 

15. 

0 . 

.7304 

1.1374 

3.2612 

3.4742 

4.4514 

6.9869 

7.8168 

15. 

5. 

.7188 

1.1413 

3.2612 

3.4858 

4.3970 

6.9762 

7.8331 

15. 

10. 

.7095 

1.1449 

3.2612 

3.4971 

4.3546 

6.9761 

7.8492 

15. 

15. 

.7023 

1.1481 

3.2611 

3.5078 

4.3221 

6.9868 

7.8650 

15. 

20. 

.6969 

1.1510 

3.2610 

3.5179 

4.2978 

7.0084 

7.8804 

15. 

25. 

.6932 

1.1535 

3.2609 

3.5273 

4.2805 

7.0408 

7.8951 

15. 

30. 

.6908 

1.1557 

3.2607 

3.5359 

4.2689 

7.0832 

7.9091 

20. 

-30. 

.8727 

1.1038 

3.2588 

3.4014 

5.3377 

7.2807 

7.7218 

20. 

-25. 

.8387 

1.1123 

3.2590 

3.4173 

5.0810 

7.2253 

7.7416 

20. 

-20. 

.8093 

1.1201 

3.2593 

3.4334 

4.8846 

7.1729 

7.7621 

20. 

-15. 

.7841 

1.1275 

3.2595 

3.4498 

4.7330 

7.1251 

7.7834 

20. 

-10. 

.7628 

1.1343 

3.2596 

3.4663 

4.6151 

7.0832 

7.8054 

20. 

-5. 

.7448 

1.1405 

3.2597 

3.4829 

4.5227 

7.0488 

7.8281 

20. 

0 . 

.7299 

1.1463 

3.2597 

3.4995 

4.4501 

7.0230 

7.8512 

20. 

5. 

.7178 

1.1515 

3.2597 

3.5159 

4.3931 

7.0068 

7.8748 

20. 

10. 

.7082 

1.1563 

3.2596 

3.5320 

4.3485 

7.0010 

7.8987 

20. 

15. 

.7008 

1.1604 

3.2595 

3.5477 

4.3142 

7.0062 

7.9227 

20. 

20. 

.6955 

1.1641 

3.2593 

3.5629 

4.2883 

7.0226 

7.9466 

20. 

25. 

.6919 

1.1673 

3.2590 

3.5775 

4.2695 

7.0502 

7.9703 

20. 

30. 

.6898 

1.1699 

3.2588 

3.5912 

4.2565 

7.0884 

7.9935 
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25. 

-30. 

.8814 

1.1015 

3.2564 

3.4082 

5.3776 

7.3415 

7.7305 

25. 

-25. 

.8446 

1.1130 

3.2568 

3.4276 

5.1067 

7.2866 

7.7550 

25. 

-20. 

.8130 

1.1235 

3.2572 

3.4477 

4.9005 

7.2333 

7.7808 

25. 

-15. 

.7862 

1.1332 

3.2574 

3.4683 

4.7422 

7.1832 

7.8080 

25. 

-10. 

.7636 

1.1420 

3.2576 

3.4894 

4.6195 

7.1379 

7.8365 

25. 

-5. 

.7447 

1.1502 

3.2578 

3.5109 

4.5237 

7.0990 

7.8665 

25. 

0 . 

.7292 

1.1575 

3.2578 

3.5326 

4.4484 

7.0680 

7.8977 

25. 

5. 

.7168 

1.1641 

3.2578 

3.5545 

4.3893 

7.0463 

7.9302 

25. 

10. 

.7070 

1.1699 

3.2576 

3.5763 

4.3430 

7.0347 

7.9638 

25. 

15. 

.6997 

1.1750 

3.2574 

3.5979 

4.3072 

7.0341 

7.9984 

25. 

20. 

.6945 

1.1794 

3.2572 

3.6192 

4.2801 

7.0449 

8.0336 

25. 

25. 

.6912 

1.1830 

3.2568 

3.6399 

4.2600 

7.0673 

8.0693 

25. 

30. 

.6896 

1.1859 

3.2564 

3.6599 

4.2460 

7.1008 

8.1051 

30. 

-30. 

.8898 

1.1003 

3.2537 

3.4199 

5.4094 

7.4085 

7.7461 

30. 

-25. 

.8500 

1.1152 

3.2542 

3.4433 

5.1268 

7.3548 

7.7758 

30. 

-20. 

.8164 

1.1285 

3.2547 

3.4676 

4.9126 

7.3011 

7.8076 

30. 

-15. 

.7879 

1.1407 

3.2551 

3.4930 

4.7489 

7.2490 

7.8416 

30. 

-10. 

.7641 

1.1518 

3.2554 

3.5193 

4.6223 

7.2004 

7.8780 

30. 

-5. 

.7444 

1.1618 

3.2555 

3.5463 

4.5237 

7.1572 

7.9169 

30. 

0 . 

.7284 

1.1709 

3.2556 

3.5740 

4.4463 

7.1211 

7.9581 

30. 

5. 

.7157 

1.1788 

3.2555 

3.6022 

4.3856 

7.0937 

8.0019 

30. 

10. 

.7059 

1.1858 

3.2554 

3.6308 

4.3380 

7.0761 

8.0481 

30. 

15. 

.6988 

1.1917 

3.2551 

3.6595 

4.3011 

7.0695 

8.0965 

30. 

20. 

.6940 

1.1966 

3.2547 

3.6880 

4.2730 

7.0744 

8.1471 

30. 

25. 

.6912 

1.2006 

3.2542 

3.7163 

4.2521 

7.0912 

8.1995 

30. 

30. 

.6902 

1.2036 

3.2537 

3.7439 

4.2374 

7.1196 

8.2532 
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Table B.2: The variation of torsion coupling ($0.75) as a function of sweep angles for 
the lag and flap modes 


A 0 

deg 

A, 

deg 

1 lag 

1 flap 

torsion 

2 flap 

2 lag 3 flap 

-30. 

-30. 

.1003 

-.1169 

.9198 

.2421 

-.0570 .0011 

-30. 

-25. 

.1046 

-.1034 

.9198 

.2273 

-.0692 .0017 

-30. 

-20. 

.1087 

-.0868 

.9199 

.2034 

-.0795 .0019 

-30. 

-15. 

.1123 

-.0674 

.9199 

.1695 

-.0873 .0018 

O 

CO 

l 

-10. 

.1152 

-.0461 

.9200 

.1249 

-.0926 .0015 

-30. 

-5. 

.1170 

-.0233 

.9200 

.0686 

-.0958 .0008 

-30. 

0 . 

.1176 

.0000 

.9200 

.0000 

-.0978 .0000 

-30. 

5. 

.1168 

.0232 

.9200 

-.0817 

-.0990 -.0010 

-30. 

10. 

.1144 

.0456 

.9200 

-.1773 

-.1000 -.0021 

-30. 

15. 

.1104 

.0664 

.9199 

-.2874 

-.1011 -.0032 

-30. 

20. 

.1045 

.0851 

.9199 

-.4125 

-.1021 -.0042 

-30. 

25. 

.0968 

.1010 

.9198 

-.5534 

-.1026 -.0052 

-30. 

30. 

.0873 

.1137 

.9198 

-.7108 

-.1024 -.0061 

-25. 

-30. 

.0906 

-.1200 

.9201 

.3080 

-.0459 .0032 

-25. 

-25. 

.0936 

-.1062 

.9201 

.2869 

-.0573 .0032 

-25. 

-20. 

.0967 

-.0890 

.9202 

.2547 

-.0673 .0030 

-25. 

-15. 

.0995 

-.0692 

.9202 

.2106 

-.0752 .0025 

-25. 

-10. 

.1019 

-.0472 

.9202 

.1539 

-.0807 .0019 

-25. 

-5. 

.1034 

-.0239 

.9203 

.0839 

-.0843 .0010 

-25. 

0 . 

.1039 

.0000 

.9203 

.0000 

-.0864 .0000 

-25. 

5. 

.1032 

.0238 

.9203 

-.0982 

-.0879 -.0011 

-25. 

10. 

.1010 

.0468 

.9202 

-.2110 

-.0894 -.0023 

-25. 

15. 

.0973 

.0682 

.9202 

-.3388 

-.0913 -.0034 

-25. 

20. 

.0920 

.0874 

.9202 

-.4817 

-.0935 -.0046 

-25. 

25. 

.0849 

.1038 

.9201 

-.6399 

-.0957 -.0056 

-25. 

30. 

.0762 

.1170 

.9201 

-.8134 

-.0975 -.0064 
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- 20 . 

- 30 . 

.0783 -.1225 

.9204 

.3873 

-.0342 

.0046 

- 20 . 

- 25 . 

.0798 -.1083 

.9204 

.3576 

-.0445 

.0043 

- 20 . 

- 20 . 

.0818 -.0908 

.9204 

.3145 

-.0538 

.0038 

- 20 . 

- 15 . 

.0838 -.0706 

.9205 

.2576 

-.0614 

.0031 

- 20 . 

- 10 . 

.0855 -.0482 

.9205 

.1864 

-.0668 

.0022 

- 20 . 

- 5 . 

.0867 -.0244 

.9205 

.1005 

-.0703 

.0011 

- 20 . 

0. 

.0871 

.0000 

.9205 

.0000 

-.0725 

.0000 

- 20 . 

5 . 

.0864 

.0243 

.9205 

-.1152 

-.0741 -.0012 

- 20 . 

10 . 

.0846 

.0479 

.9205 

-.2449 

-.0759 -.0024 

- 20 . 

15 . 

.0813 

.0698 

.9205 

-.3886 

-.0784 -.0037 

- 20 . 

20 . 

.0766 

.0895 

.9204 

-.5456 

-.0816 -.0048 

- 20 . 

25 . 

.0704 

.1064 

.9204 

-.7151 

-.0855 -.0058 

- 20 . 

30 . 

.0628 

.1199 

.9204 

-.8959 

-.0893 -.0067 

- 15 . 

- 30 . 

.0638 -.1244 

.9206 

.4802 

-.0223 

.0056 

- 15 . 

- 25 . 

.0638 -.1100 

.9206 

.4386 

-.0312 

.0050 

- 15 . 

- 20 . 

.0644 -.0922 

.9206 

.3816 

-.0394 

.0043 

- 15 . 

- 15 . 

.0655 -.0717 

.9206 

.3090 

-.0462 

.0034 

- 15 . 

- 10 . 

.0666 -.0490 

.9207 

.2209 

-.0512 

.0024 

- 15 . 

- 5 . 

.0674 -.0248 

.9207 

.1177 

-.0544 

.0012 

- 15 . 

0. 

.0677 

.0000 

.9207 

.0000 

-.0564 

.0000 

- 15 . 

5 . 

.0671 

.0248 

.9207 

-.1314 

-.0579 -.0013 

- 15 . 

10 . 

.0656 

.0487 

.9207 

-.2753 

-.0598 -.0026 

- 15 . 

15 . 

.0629 

.0711 

.9206 

-.4302 

-.0627 -.0038 

- 15 . 

20 . 

.0590 

.0913 

.9206 

-.5945 

-.0669 -.0050 

- 15 . 

25 . 

.0538 

.1085 

.9206 

-.7658 

-.0722 -.0060 

- 15 . 

30 . 

.0473 

.1225 

.9206 

-.9417 

-.0780 -.0069 

- 10 . 

- 30 . 

.0477 -.1258 

.9208 

.5848 

-.0106 

.0062 

- 10 . 

- 25 . 

.0459 -.1113 

.9208 

.5272 

-.0176 

.0055 

- 10 . 

- 20 . 

.0452 -.0933 

.9208 

.4525 

-.0244 

.0047 

- 10 . 

- 15 . 

.0452 -.0726 

.9208 

.3613 

-.0301 

.0037 

- 10 . 

- 10 . 

.0457 -.0496 

.9208 

.2546 

-.0343 

.0025 

- 10 . 

- 5 . 

.0461 -.0252 

.9208 

.1336 

-.0370 

.0013 

- 10 . 

0. 

.0462 

'.0000 

.9208 

.0000 

-.0385 

.0000 

- 10 . 

5 . 

.0459 

.0251 

.9208 

-.1444 

-.0398 -.0013 

- 10 . 

10 . 

.0447 

.0494 

.9208 

-.2974 

-.0417 -.0027 

- 10 . 

15 . 

.0427 

.0722 

.9208 

-.4565 

-.0448 -.0039 

o 

r-H 

1 

20 . 

.0396 

.0927 

.9208 

-.6188 

-.0497 -.0051 

- 10 . 

25 . 

.0354 

.1103 

.9208 

-.7810 

-.0562 -.0061 

- 10 . 

30 . 

.0303 

.1245 

.9208 

-.9395 

-.0639 -.0070 
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- 5 . 

- 30 . 

.0303 

-.1266 

.9209 

.6958 

.0006 

.0066 

- 5 . 

- 25 . 

.0267 

-.1120 

.9209 

.6171 

-.0043 

.0059 

- 5 . 

- 20 . 

.0246 

-.0940 

.9209 

.5209 

-.0093 

.0049 

- 5 . 

- 15 . 

.0237 

-.0731 

.9209 

.4089 

-.0136 

.0038 

- 5 . 

- 10 . 

.0234 

-.0500 

.9209 

.2832 

-.0168 

.0026 

- 5 . 

- 5 . 

.0234 

-.0254 

.9209 

.1460 

-.0186 

.0013 

- 5 . 

0 . 

.0235 

.0000 

.9209 

.0000 

-.0196 

.0000 

- 5 . 

5 . 

.0233 

.0254 

.9209 -.1520 

-.0204 -.0014 

- 5 . 

10 . 

.0225 

.0499 

.9209 -.3070 

-.0219 -.0027 

- 5 . 

15 . 

.0212 

.0729 

.9209 -.4618 

-.0251 -.0040 

- 5 . 

20 . 

.0190 

.0937 

.9209 -.6129 

-.0305 -.0051 

- 5 . 

25 . 

.0160 

.1116 

.9209 -.7566 

-.0381 -.0062 

- 5 . 

30 . 

.0123 

.1260 

.9209 -.8893 

-.0475 -.0070 

0 . 

- 30 . 

.0123 

-.1267 

.9209 

.8025 

.0111 

.0069 

0 . 

- 25 . 

.0068 

-.1122 

.9209 

.6980 

.0085 

.0061 

0 . 

- 20 . 

.0033 

-.0941 

.9209 

.5778 

.0055 

.0051 

0 . 

- 15 . 

.0014 

-.0732 

.9209 

.4447 

.0028 

.0039 

0 . 

- 10 . 

.0004 

-.0501 

.9209 

.3019 

.0010 

.0027 

0 . 

- 5 . 

.0000 

-.0254 

.9209 

.1526 

.0001 

.0014 

0 . 

0 . 

.0000 

.0000 

.9209 

.0000 

.0000 

.0000 

0 . 

5 . 

.0000 

.0254 

.9209 -.1526 

-.0002 -.0014 

0 . 

10 . 

-.0003 

.0501 

.9209 -.3019 

-.0013 -.0027 

0 . 

15 . 

-.0010 

.0732 

.9209 -.4447 

-.0043 -.0039 

0 . 

20 . 

-.0022 

.0941 

.9209 -.5778 

-.0099 -.0051 

0 . 

25 . 

-.0039 

.1122 

.9209 -.6980 

-.0183 -.0061 

0 . 

30 . 

-.0062 

.1267 

.9209 -.8025 

-.0292 -.0069 

5 . 

- 30 . 

-.0059 

-.1259 

.9209 

.8893 

.0206 

.0070 

5 . 

- 25 . 

-.0132 

-.1115 

.9209 

.7566 

.0204 

.0062 

5 . 

- 20 . 

-.0180 

-.0936 

.9209 

.6128 

.0196 

.0051 

5 . 

- 15 . 

-.0210 

-.0729 

.9209 

.4618 

.0187 

.0040 

5 . 

- 10 . 

-.0226 

-.0499 

.9209 

.3070 

.0184 

.0027 

5 . 

- 5 . 

-.0233 

-.0253 

.9209 

.1520 

.0187 

.0014 

5 . 

0 . 

-.0235 

.0000 

.9209 

.0000 

.0196 

.0000 

5 . 

5 . 

-.0233 

.0254 

.9209 -.1460 

.0202 

-.0013 

5 . 

10 . 

-.0231 

.0500 

.9209 -.2832 

.0197 

-.0026 

5 . 

15 . 

-.0231 

.0731 

.9209 -.4089 

.0170 -.0038 

5 . 

20 . 

-.0233 

.0940 

.9209 -.5209 

.0114 

-.0049 

5 . 

25 . 

-.0238 

.1120 

.9209 -.6171 

.0024 -.0059 

5 . 

30 . 

-.0246 

.1266 

.9209 -.6958 

-.0097 -.0066 
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10. 

-30. 

-.0236 

-.1239 

.9208 

.9394 

.0290 

.0070 

10. 

-25. 

-.0328 

-.1098 

.9208 

.7810 

.0312 

.0061 

10. 

-20. 

-.0389 

-.0923 

.9208 

.6188 

.0326 

.0051 

10. 

-15. 

-.0427 

-.0719 

.9208 

.4565 

.0337 

.0039 

10. 

-10. 

-.0450 

-.0492 

.9208 

.2974 

.0350 

.0027 

10. 

-5. 

-.0460 

-.0250 

.9208 

.1444 

.0367 

.0013 

10. 

0 . 

-.0462 

.0000 

.9208 

.0000 

.0385 

.0000 

10. 

5. 

-.0459 

.0251 

.9208 -.1336 

.0401 -.0013 

10. 

10. 

-.0453 

.0494 

.9208 -.2546 

.0403 -.0025 

10. 

15. 

-.0445 

.0723 

.9208 -.3613 

.0382 -.0037 

10. 

20. 

-.0437 

.0930 

.9208 -.4525 

.0327 -.0047 

10. 

25. 

-0430 

.1109 

.9208 -.5273 

.0235 -.0055 

10. 

30. 

-.0423 

.1255 

.9208 -.5849 

.0104 -.0062 

15. 

-30. 

-.0405 

-.1203 

.9206 

.9415 

.0361 

.0069 

15. 

-25. 

-.0513 

-.1068 

.9206 

.7658 

.0408 

.0060 

15. 

-20. 

-.0585 

-.0899 

.9206 

.5945 

.0444 

.0050 

15. 

-15. 

-.0632 

-.0701 

.9206 

.4303 

.0475 

.0038 

15. 

-10. 

-.0660 

-.0481 

.9207 

.2753 

.0504 

.0026 

15. 

-5. 

-.0674 

-.0245 

.9207 

.1314 

.0534 

.0013 

15. 

0 . 

-.0677 

.0000 

.9207 

.0000 

.0564 

.0000 

15. 

5. 

-0672 

.0245 

.9207 -.1177 

.0588 -.0012 

15. 

10. 

-.0661 

.0484 

.9207 -.2209 

.0599 -.0024 

15. 

15. 

-.0646 

.0709 

.9206 -.3090 

.0584 -.0034 

15. 

20. 

-.0629 

.0913 

.9206 -.3817 

.0534 -.0043 

15. 

25. 

-.0610 

.1089 

.9206 -.4387 

.0441 -.0050 

15. 

30. 

-.0589 

.1233 

.9206 -.4803 

.0304 -.0056 

20. 

-30. 

-.0559 

-.1150 

.9204 

.8956 

.0420 

.0067 

20. 

-25. 

-.0682 

-.1025 

.9204 

.7149 

.0489 

.0058 

20. 

-20. 

-.0764 

-.0864 

.9204 

.5456 

.0547 

.0048 

20. 

-15. 

-.0818 

-.0675 

.9205 

.3886 

.0597 

.0037 

20. 

-10. 

-.0851 

-.0464 

.9205 

.2450 

.0642 

.0025 

20. 

-5. 

-.0867 

-.0236 

.9205 

.1153 

.0685 

.0012 

20. 

0 . 

-.0871 

.0000 

.9205 

.0000 

.0725 

.0000 

20. 

5. 

-.0864 

.0238 

.9205 -.1006 

.0759 -.0011 

20. 

10. 

-.0849 

.0470 

.9205 -.1865 

.0778 -.0022 

20. 

15. 

-.0828 

.0688 

.9205 -.2578 

.0771 -.0031 

20. 

20. 

-.0802 

.0887 

.9204 -.3147 

.0728 -.0038 

20. 

25. 

-.0772 

.1060 

.9204 -.3578 

.0638 -.0043 

20. 

30. 

-.0738 

.1200 

.9204 -.3876 

.0497 -.0046 
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25. 

-30. 

-.0695 

-.1079 

.9201 

.8130 

.0466 

.0064 

25. 

-25. 

-.0831 

-.0967 

.9201 

.6398 

.0555 

.0056 

25. 

-20. 

-.0920 

-.0819 

.9202 

.4818 

.0632 

.0046 

25. 

-15. 

-.0979 

-.0642 

.9202 

.3389 

.0700 

.0035 

25. 

-10. 

-.1016 

-.0442 

.9202 

.2112 

.0760 

.0023 

25. 

-5. 

-.1035 

-.0226 

.9203 

.0983 

.0814 

.0011 

25. 

0 . 

-.1039 

.0000 

.9203 

.0000 

.0864 

.0000 

25. 

5. 

-.1031 

.0228 

.9203 

-.0840 

.0907 

-.0010 

25. 

10. 

-.1012 

.0451 

.9202 

-.1541 

.0934 

-.0019 

25. 

15. 

-.0985 

.0661 

.9202 

-.2109 

.0938 

-.0025 

25. 

20. 

-.0951 

.0853 

.9202 

-.2551 

.0903 

-.0030 

25. 

25. 

-.0911 

.1020 

.9201 

-.2874 

.0819 

-.0032 

25. 

30. 

-.0865 

.1157 

.9201 

-.3086 

.0679 

-.0032 

30. 

-30. 

-.0809 

-.0988 

.9198 

.7103 

.0499 

.0061 

30. 

-25. 

-.0953 

-.0895 

.9198 

.5534 

.0605 

.0052 

30. 

-20. 

-.1048 

-.0764 

.9199 

.4127 

.0699 

.0043 

30. 

-15. 

-.1111 

-.0601 

.9199 

.2877 

.0782 

.0032 

30. 

-10. 

-.1151 

-.0415 

.9200 

.1776 

.0855 

.0021 

30. 

-5. 

-.1172 

-.0213 

.9200 

.0819 

.0920 

.0010 

30. 

0 . 

-.1176 

.0000 

.9200 

.0000 

.0978 

.0000 

30. 

5. 

-.1167 

.0216 

.9200 

-.0688 

.1027 

-.0008 

30. 

10. 

-.1145 

.0427 

.9200 

-.1252 

.1064 

-.0015 

30. 

15. 

-.1113 

.0628 

.9199 

-.1701 

.1078 

-.0018 

30. 

20. 

-.1072 

.0812 

.9199 

-.2040 

.1056 

-.0019 

30. 

25. 

-.1023 

.0972 

.9198 

-.2281 

.0983 

-.0017 

30. 

30. 

-.0968 

.1104 

.9198 

-.2430 

.0848 

-.0011 





